AD-A115  703  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/G  15/7 

AN  INVESTIGATION  OF  THE  BOMBER  ANO  TANKER  MATING  PROCESS  IN  THE— ETC(U) 
MAR  62  ML  MACELHANEY.  J  W  STANFIELD 

UNCLASSIFIED  AFI T/G5T/QS/62M-B _ NL _ 


OTIC  F![[  COPY  AD  All  5*7  03 


APIT/GST/OS/82M-8 


AN  INVESTIGATION  OP  THE  BOMBER  AND 
TANKER  MATING  PROCESS  IN  THE 
SINGLE  INTEGRATED  OPERATIONS  PLAN 


THESIS 

William  L.  MacElhaney,  Capt,  USAF 
James  W.  Stanfield,  Capt,  USAF 

AF IT/GST/OS/8  2M-8 


DTIC 


ELECTE 
JUN  18  1982 


D 


Approved  for  public  release;  distribution  unlimited 


"Original  contains  color  ", 
plates :  All  DTIC  reproduct, 
ions  will  be  in  black  and 

i t  e  * 

82  06  16  079 


AFIT/ GST/OS/8  2M-8 


"Original  contains  oolor 
plates:  All  D1IC  reproduct¬ 
ions  will  be  in  black  and 
uh i t  e • 


Accession  Fop 

NT IS  GRA4I 
DTIC  TAB 
Unannounced 
Justification- 


X" 

n 

n 


By  1 

__Distr 

Aval 

Dist 

ft 

ibutioj 

labilll 

Avail 

Spec 

a/ 

-y  Codes 

and/or 

ial 

AN  INVESTIGATION  OF  THE  BOMBER  AND  TANKER  MATING 
PROCESS  IN  THE  SINGLE  INTEGRATED  OPERATIONS  PLAN 


THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


by 


William  L.  MacElhaney,  BS  James  W.  Stanfield,  BS 

Capt  USAF  Capt  USAF 

Graduate  Strategic  and  Tactical  Sciences 
March  1982 


Approved  for  public  release;  distribution  unlimited 


Preface 


The  process  of  assigning  tankers  to  bombers  in  the 
Single  Integrated  Operations  Plan  is  very  critical.  Tanker 
shortages  require  that  these  resources  be  utilized  in  the 
best  way  possible.  The  goal  of  this  thesis  effort  was  to 
find  a  way  of  improving  on  the  current  methods  used  in 
assigning  tankers,  and  hopefully  we  have  done  so. 

This  goal  could  not  have  been  attained  without  the 
help  of  many  people,  and  we  wish  to  thank  all  of  those  that 
helped  us  in  any  way. 

In  particular,  we  thank  Capt  Jeff  Goodlett,  HQ  SAC/ 
XOXF,  for  his  enthusiastic  help  in  providing  us  with  all 
of  the  required  data  and  planning  documents  used  in  our 
research  of  the  problem.  We  thank  Dr.  Ken  Kast  of  the 
Logicon  Corporation  for  providing  us  with  the  methodology 
employed  in  the  current  Mating  and  Ranging  Program,  as  well 
as  for  giving  us  his  thoughts  on  our  approach  to  the 
problem.  We  thank  our  thesis  advisor.  Major  Gerald  R. 
Armstrong,  and  our  reader,  Major  Ivy  D.  Cook,  for  their 
helpful  suggestions  and  conscientious  guidance  throughout 
the  preparation  of  this  thesis. 

Finally,  and  most  of  all,  we  wish  to  thank  our 
families  for  providing  us  with  the  encouragement  and  sup¬ 
port  needed  to  complete  a  long  and  exhausting  effort. 
Special  thanks  goes  to  Barbara  Stanfield  for  her  many 
hours  spent  in  the  typing  of  rough  drafts. 
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Abstract 

The  survivability  of  the  strategic  bomber  force 
during  Emergency  War  Order  missions  is  of  primary  concern 
to  the  Strategic  Air  Command.  Since  the  Soviet  Union  and 
its  Warsaw  Pact  allies  possess  the  most  extensive  air 
defense  networks  in  the  world,  a  penetrating  bomber  force 
must  fly  as  low  as  possible  for  as  long  as  possible.  This 
tactic  minimizes  probability  of  detection  and  vulnerability 
to  defensive  threats.  It  also  greatly  increases  the  fuel 
required  to  complete  the  mission.  This  additional  fuel  is 
supplied  by  one  or  more  in-flight  refuelings. 

'"•The  initial  objective  of  this  thesis  was  to  develop 
a  method  for  assigning  tankers  to  the . bomber  force  in  an 
optimal  manner.  As  the  study  progressed  however,  it  became 
clear  that  obtaining  a  truly  optimal  solution  using  mathe¬ 
matical  programming  techniques  cannot  be  guaranteed  due  to 
the  nature  and  complexity  of  the  problem.  As  a  result  the 
emphasis  of  the  study  was  shifted  to  developing  an  improved 
method  for  solving  the  problem. 

Two  heuristic  methods  were  investigated.  The  first 
method  used  network  theory  in  an  attempt  to  minimize  the 
costs  of  assigning  tankers  to  the  bombers.  The  second 
method  was  based  on  the  so-called  "greedy"  method.  This 


viii 


method  basically  made  the  assignments  in  the  order  of 
decreasing  marginal  cost  improvements.  These  two  methods 
were  evaluated  against  each  other  and  the  current  method 
by  means  of  several  example  problems.  Both  methods  yielded 
better  results  than  the  one  currently  in  use,  with  the  net¬ 
work  method  appearing  to  be  the  best. 


AN  INVESTIGATION  OF  THE  BOMBER  AND  TANKER 
MATING  PROCESS  IN  THE  SINGLE 
INTEGRATED  OPERATIONS  PLAN 


I.  Introduction 


Background 

The  military  forces  of  the  United  States  seek  to 
deter  aggression  from  other  nations  by  maintaining  forces 
capable  of  responding  to  threats  across  the  spectrum  of  war¬ 
fare.  The  United  States'  strategic  nuclear  forces  deter 
war  in  general  and  nuclear  war  in  particular,  by  maintain¬ 
ing  a  TRIAD  of  forces  that  include  land-based  interconti¬ 
nental  ballistic  missiles  (ICBM),  submarine-launched  ballis¬ 
tic  missiles  (SLBM) ,  and  manned  bombers.  V 

The  manned  bomber  leg  of  the  TRIAD  consists  of  B-52 
and  FB-111  bomber  aircraft  which  are  supported  by  KC-135 
tanker  aircraft.  The  tanker  aircraft  add  range,  flexi¬ 
bility,  and  responsiveness  to  the  bomber  forces  and  allow 
them  to  fly  long-range  strike  missions  deep  into  enemy 
territory  and  return  home  or  to  post-strike  bases  in 
friendly  territory  (Ref  11:45) . 

The  general  profile  of  the  bomber  missions  is  to 
takeoff,  fly  at  high  altitude  (to  include  air  refueling) 
until  reaching  the  perimeter  of  enemy  defensive  coverage. 


descend  to  low  altitude  to  avoid  defenses ,  and  strike 
assigned  targets.  After  striking  all  targets,  the  bomber 
exits  enemy  defensive  coverage  at  low  altitude  and  then 
climbs  back  to  high  altitude  for  the  flight  home  or  to 
post-strike  bases  in  friendly  territory. 

In  order  to  accomplish  the  entire  mission  as  out¬ 
lined  above,  a  bomber  may  require  one,  two,  or  even  three 
air  refuelings  prior  to  descending  to  low  level.  These 
refuelings  are  provided  by  tanker  aircraft  that  may  be 
co-located  with  the  bomber  or  located  at  another  base.  In 
the  latter  case,  the  bomber  and  tanker  rendezvous  along 
the  bomber's  route  of  flight.  After  offloading  fuel  to  the 
bomber,  the  tanker  recovers  at  a  pre-determined  post¬ 
refueling  base  (PRB) .  This  complete  process  is  illustrated 
in  Figure  1. 

The  bomber's  air  refueling  requirements  are  deter¬ 
mined  by  flying  its  mission  in  reverse.  This  is  accom¬ 
plished  by  starting  at  the  post-strike  base  and  working 
backwards  to  a  point  just  prior  to  where  the  bomber  has  to 
descend  to  low  altitude.  This  point  is  designated  as  the 
entry  point  (EP) .  The  difference  between  the  bomber's 
fuel  at  the  EP  without  any  air  refueling  and  the  fuel 
actually  needed  to  complete  the  mission  is  then  used  to 
determine  the  onload  and  the  number  of  tankers  required  to 
deliver  it.  Any  onload  less  than  this  figure  will  require 
the  bomber  to  use  degraded  tactics,  i.e.,  descend  to  low 
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Departure  Base 


Post  Strike  Base 


Bomber  Altitude  Profile 


Departure  Base  PRB 

Tanker  Altitude  Profile 


Tanker-Bomber  Route  Profile 
Figure  1.  Bomber/Tanker  Route  Structure 


altitude  later  and/or  climb  out  of  low-level  earlier  than 
desired.  Unfortunately,  this  is  the  case  for  a  large  pro¬ 
portion  of  the  bomber  force  because  there  are  not  enough 
tanker  aircraft  to  provide  all  of  the  required  fuel  (Ref  8) 

Current  Procedures 

Since  tanker  resources  are  limited  and  bomber  fuel 
requirements  are  so  high,  it  is  extremely  important  to 
utilize  the  available  tankers  in  the  most  efficient  manner; 
however,  this  is  easier  said  than  done.  Under  full  scale 
implementation  of  the  Single  Integrated  Operations  Plan 
(SIOP) ,  hundreds  of  bombers  require  refuelings  from  an  even 
larger  number  of  tankers  that  are  deployed  in  numerous  loca 
tions,  and  subsequently  recover  at  different  locations. 

The  fuel  available  from  each  of  these  tankers  depends  on 
where  it  comes  from,  where  the  refueling  occurs,  and  where 
it  recovers  after  refueling. 

Prior  to  1980,  the  mission  planners  at  SAC  Head¬ 
quarters  used  manual  procedures  to  accomplish  bomber/tanker 
mating.  This  process  involved  continuous  iterations  until 
a  feasible  solution  was  obtained.  It  was  a  very  arduous 
and  time-consuming  procedure,  sometimes  requiring  several 
weeks  to  complete  (Ref  8) .  Due  to  the  complexity  of  the 
problem  and  the  emphasis  on  a  feasible  solution,  there  was 
no  assurance  that  the  solution  obtained  was  even  close  to 
optimum. 


Starting  in  1980 ,  the  mission  planners  began  to  use 
a  computerized  algorithm  developed  by  the  Logicon  Corpora¬ 
tion  to  assist  them  in  the  mating  process  (Ref  8) .  This 
Mating  and  Ranging  Program  (MARP)  performs  the  entire 
mission-planning  process  including  assigning  tankers  to 
bombers  and  determining  the  air  refueling  locations.  The 
assignments  are  determined  by  flowing  the  tankers  through 
a  network,  the  details  of  which  will  be  outlined  in 
Chapter  III.  This  program  has  aided  the  mission  planners 
a  great  deal,  but  some  manual  calculations  and  matings  are 
still  required.  Further  efforts  at  improving  the  assign¬ 
ment  process  have  met  with  little  success  to  date  (Ref  9)  . 

Problem  Statement 

As  a  result  of  tanker  shortages,  a  large  number  of 
bomber  sorties  must  resort  to  degraded  tactics  in  order  to 
reach  their  post-strike  bases  with  required  fuel  reserves. 
The  introduction  of  air  launched  cruise  missiles  will  fur¬ 
ther  increase  bomber  fuel  requirements.  The  problem  then 
is  to  utilize  the  tanker  force  in  the  best  possible  manner. 
Ideally,  the  tankers  should  be  assigned  so  as  to  meet  all 
of  the  bomber  EP  fuel  requirements.  Since  the  tanker 
shortage  precludes  this  possibility,  they  should  then  be 
assigned  to  meet  as  many  of  the  requirements  as  possible 
while  minimizing  the  shortages  of  those  bombers  whose 
requirements  can  not  be  met. 
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Objective 

The  objective  of  this  research  effort  was  to 
develop  an  optimal  or  near  optimal  methodology  for  mating 
bombers  and  tankers  in  the  SIOP.  The  goals  of  this  method- 
ology  are  to  reduce  the  number  of  bombers  requiring  degraded 
tactics  and/or  reduce  the  duration  of  these  tactics. 

Two  different  methods  were  investigated  for  accom¬ 
plishing  this  objective.  The  first  method  uses  a  network 
algorithm  to  minimize  the  costs  associated  with  assigning 
tankers  to  the  bombers.  This  method  was  developed  indepen¬ 
dently  of  the  Logicon  method  and  was  believed  to  hold  the 
most  promise  for  obtaining  an  optimal  solution.  The  second 
method  is  based  on  the  "greedy"  method.  This  method  assigns 
tankers  based  on  their  marginal  cost  contributions.  This 
method  was  used  because  it  is  easily  implemented  and  has 
been  applied  to  a  wide  variety  of  problems  (Ref  17:59-70). 

It  therefore  served  the  additional  purpose  of  being  a  bench¬ 
mark  against  which  both  the  Logicon  and  network  methods 
could  be  evaluated. 

Scope  and  Limitations 

As  previously  mentioned,  the  full  SIOP  involves 
hundreds  of  bombers  and  tankers  and  numerous  bases.  In 
order  to  reduce  the  problem  to  manageable  proportions,  this 
study  was  confined  to  normal  day-to-day  alert  force  air¬ 
craft.  Limiting  the  problem  in  this  manner  reduces  the 


number  of  aircraft  involved  to  under  300.  Any  methodology 
developed  to  handle  this  problem  could  subsequently  be 
expanded  to  deal  with  the  larger  problem. 

Aircraft  included  in  the  study  are  limited  to  the 
B-52H  and  the  KC-135A.  This  limitation  was  imposed  to  take 
advantage  of  the  authors'  knowledge  of  these  specific  air¬ 
craft  and  to  eliminate  the  complexities  involved  with 
tracking  five  different  types  of  aircraft  (B-52D,  B-52G, 
B-52H,  FB-111,  and  KC-135A)  and  their  different  performance 
computations.  This  limitation  should  not  affect  the  overall 
solution  because  the  basic  requirements  remain  the  same 
regardless  of  the  aircraft  type. 

In  addition  to  mating  bombers  and  tankers,  the 
Logicon  program  performs  numerous  other  functions  such  as 
conflict  resolutions,  avoiding  flight  over  major  target 
complexes,  and  providing  detailed  flight  plans  (Ref  8) . 

These  functions  were  beyond  the  scope  of  this  effort  due 
to  time  and  manpower  constraints  and  were  not  considered. 

Finally,  most  of  the  data  involving  SIOP  forces 
is  classified  at  the  SECRET  level  or  higher.  It  is  for 
this  reason  that  most  figures  are  quoted  as  approximations 
only.  Furthermore,  the  methodologies  are  developed  and 
evaluated  using  a  combination  of  hypothetical  and  real 
numbers  and  locations  to  avoid  the  problems  of  classifica¬ 
tion.  At  the  same  time,  every  effort  has  been  made  to 
make  the  various  example  problems  as  realistic  as  possible. 


Sequence  of  Presentation 

The  remainder  of  this  thesis  is  devoted  to  the 
accomplishment  of  the  objective  which  was  stated  earlier. 
Chapter  II  details  the  prototype  problem  used  to  develop 
the  methodologies  and  computer  programs  for  the  proposed 
new  methods.  It  also  outlines  the  assumptions  used  in 
developing  and  evaluating  them.  The  theory  and  methodology 
of  each  of  the  new  methods  as  well  as  the  current  method 
(MARP)  are  discussed  in  Chapter  III.  Each  method  is  then 
used  to  solve  the  prototype  problem  and  four  additional 
mating  problems.  The  results  are  reported  and  analyzed 
in  Chapter  IV.  Finally,  conclusions  and  recommendations 
are  presented  in  Chapter  V. 


II.  The  Prototype  Problem 


Description 

As  noted  earlier,  this  study  was  limited  to  the 
problem  of  mating  bomber  and  tanker  aircraft  that  are  on 
normal,  day-to-day  alert.  This  limitation  reduced  the 
number  of  aircraft  involved  considerably,  but  it  did  little 
to  diminish  the  complexity  of  the  problem.  Several  hundred 
bombers  and  tankers  from  numerous  locations  still  have  to 
be  mated  for  air  refueling,  followed  by  recovery  of  the 
tankers  at  various  different  locations. 

Investigating  a  problem  of  this  magnitude  from  the 

start  would  have  been  extremely  difficult.  It  is  for  this 

reason  that  the  approach  to  modeling  advocated  by  William  T. 

Morris  was  followed  (Ref  16:707).  He  states: 

The  process  of  model  development  may  be  usually  viewed 
as  a  process  of  enrichment  or  elaboration.  One  begins 
with  very  simple  models,  quite  distinct  from  reality 
and  attempts  to  move  in  evolutionary  fashion  toward 
more  elaborate  models  which  more  nearly  reflect  the 
complexity  of  the  actual  management  situation. 

This  was  the  process  followed  in  the  conceptualization, 
modeling,  and  computerization  of  this  study.  The  initial 
model  or  prototype  problem  consisted  of  four  bombers, 
seven  tankers,  and  three  recovery  bases.  In  addition,  the 
takeoff  gross  weights  of  both  aircraft  and  their  respec¬ 
tive  fuel  consumption  rates  were  assumed  to  be  constant. 
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The  first  "enrichment  and  elaboration"  of  the  prototype 
problem  deleted  all  but  one  of  the  assumptions.  Only  the 
bomber's  takeoff  gross  weight  remained  fixed.  The  tanker's 
takeoff  gross  weight  and  the  fuel  consumption  rates  of  both 
aircraft  were  allowed  to  vary.  This  version  of  the  problem 
was  then  further  "enriched  and  elaborated"  until  the  final 
problem  consisted  of  90  bombers,  135  tankers,  and  18  post¬ 
refueling  bases. 

The  prototype  problem  is  illustrated  in  Figure  2. 
There  is  one  bomber  and  tanker  at  each  of  the  first  three 
bomber  bases  and  one  bomber  only  at  the  fourth  bomber  base. 
The  tanker-only  bases  have  two  tankers  assigned  to  them. 

The  bombers  fly  the  routes  as  depicted  from  their  bases  to 
the  entry  points.  Enroute  they  receive  one  or  more  air 
refuelings  which  are  conducted  on  the  route  segments 
denoted  by  the  small  squares.  In  this  particular  scenario, 
bombers  one  through  four  receive  one,  three,  one,  and  two 
refuelings  respectively.  After  refueling  is  completed, 
the  tankers  land  at  the  appropriate  PRB .  PRBs  one  through 
three  have  capacities  of  two,  two,  and  ten  tankers  respec¬ 
tively.  These  capacities  arise  from  ramp  space  limita¬ 
tions  and/or  servicing  capabilities. 

The  problem  then  becomes  how  to  assign  the  tankers 
to  the  refueling  tracks  and  thence  to  the  recovery  bases 
so  as  to  minimize  the  number  of  bombers  that  are  short  of 
their  entry  point  fuel  requirement  and  to  minimize  any 
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AR  -  Air  Refueling 
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TB  -  Tanker  Base 

Prototype  Problem 


shortages.  The  three  methods  discussed  in  Chapter  I 
attempt  to  do  just  this.  Their  success  in  doing  so  is 
evaluated  against  the  prototype  problem  just  described  as 
well  as  against  four  additional  problems.  Each  of  the  sue 
ceeding  problems  or  models  is  larger  than  its  predecessor , 
culminating  with  the  one  that  consists  of  90  bombers,  135 
tankers,  and  18  recovery  bases. 

Flight  Planning 

Each  of  the  methods  under  investigation  in  this 
study  attempt  to  solve  the  problem  outlined  above  by 
maximizing  the  tanker  offload  capabilities  in  one  manner 
or  another.  The  tanker  offload  capabilities,  in  turn,  are 
maximized  by  minimizing  the  fuel  required  to  deliver  them. 
The  fuel  requirements  for  delivering  these  offloads  are 
in  effect,  the  costs  of  this  study.  To  minimize  these 
costs  it  is  first  necessary  to  compute  them.  Computing 
these  figures  also  requires  additional  computations  to 
determine  fuel  consumption  rates  and  onload  and  offload 
capabilities.  The  process  of  computing  these  various  fuel 
figures  is  known  as  flight  planning.  This  process  will  be 
discussed  in  detail  in  subsequent  sections.  Prior  to  this 
however,  it  is  necessary  to  review  some  of  the  underlying 
assumptions  that  apply  to  the  mating  problem  and  the 
overall  flight  planning  process.  Any  additional 


assumptions  that  apply  to  a  particular  phase  of  the  flight 
planning  process  are  listed  under  that  phase. 


1.  The  number  of  bombers  and  tankers  and  their 
respective  bases  are  fixed. 

2.  The  available  tanker  recovery  bases  and  their 
ramp  capacities  are  fixed. 

3.  The  entry  point  and  the  fuel  required  at  that 
point  are  fixed. 

4.  All  aircraft  launch  or  takeoff  from  a  ground 
alert  posture. 

5.  All  aircraft  takeoff  at  the  same  time  under  an 
attack  warning. 

6.  Since  the  bomber  aircraft  are  not  performance 
limited,  they  takeoff  at  their  maximum  allowable  gross 
weights  in  accordance  with  standard  Emergency  War  Order 
(EWO)  planning  factors . 

7.  The  tanker  aircraft  on  the  other  hand,  are 
frequently  performance  limited  because  of  field  and  cli¬ 
matic  conditions.  Therefore  some  aircraft  takeoff  at  less 
than  maximum  allowable  gross  weight.  This  is  also  in 
accordance  with  standard  EWO  planning  factors. 

8.  All  aircraft  fly  great  circle  routes  (most 
direct)  from  their  departure  bases  to  their  entry  points. 
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refueling  tracks,  and  recovery  bases  as  applicable.  Avoid¬ 
ance  of  target  complexes  and  possible  route  conflicts  is 
not  considered. 

9.  All  flight  planning  computations  are  based  on 
standard  EWO  planning  factors  or  the  respective  aircraft 
performance  manuals  as  applicable  (Refs  4,  13). 

10.  Standard  day  conditions  were  assumed  to  apply 
throughout.  The  primary  factor  involved  in  this  assumption 
is  temperature .  Temperatures  warmer  than  standard  generally 
reduce  aircraft  performance  while  colder  temperatures 
usually  enhance  performance .  This  assumption  should  not 
affect  the  overall  results  because  it  applies  equally  to 
all  three  methods. 

11.  All  flight  planning  calculations  are  based  on 
no-wind  conditions.  Headwinds  adversely  affect  range  and 
timing  considerations  while  tailwinds  enhance  them.  This 
assumption,  like  standard  day  conditions,  should  not  affect 
the  final  results  since  it  also  applies  equally  to  all 
three  methods. 

12.  The  first  possible  refueling  point  occurs 
after  both  the  bomber  and  tanker  have  leveled  off  and 
refueling  must  be  completed  prior  to  the  entry  point. 

13.  A  bomber  will  never  delay  enroute  to  meet  a 
tanker.  This  effectively  requires  a  bomber's  tanker  to  be 
co-located  with  the  bomber  or  to  be  located  forward  of  the 
bomber's  route  of  flight.  For  example,  the  tanker  located 
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at  tanker  base  3  of  Figure  2  can  not  refuel  any  of  the 
other  bombers,  because  they  would  have  to  delay  enroute 
so  that  the  tanker  could  join  them.  On  the  other  hand, 
tankers  from  tanker  bases  1,  2,  or  4  can  refuel  bomber  2 
because  they  are  located  forward  of  the  bomber 1 s  route  of 
flight.  The  test  as  to  whether  a  particular  tanker  can 
refuel  a  given  bomber  is  based  upon  arrival  time  at  the 
start  refueling  point  which  is  designated  as  the  air  refuel¬ 
ing  control  point  (ARCP) .  If  the  tanker  cam  arrive  at  the 
ARCP  at  or  prior  to  the  bomber's  arrival  time,  the  refuel¬ 
ing  is  feasible.  If  not,  the  refueling  is  infeasible. 

14 .  If  a  tanker  arrives  at  the  ARCP  prior  to  its 
bomber,  it  will  enter  a  holding  pattern  at  the  ARCP  and 
await  the  bomber. 

15.  The  maximum  number  of  refuelings  for  a  bomber 
is  three. 

Flight  Planning  Process 

For  the  purposes  of  this  study,  the  flight  planning 
process  has  been  divided  into  four  phases.  These  phases 
are  tadceoff  to  the  ARCP,  holding  at  the  ARCP,  air  refuel¬ 
ing,  and  post-air  refueling.  These  phases  are  illustrated 
in  Figure  3  and  discussed  below. 

Takeoff  to  the  ARCP.  Takeoff,  climb,  level-off, 
and  enroute  cruise  to  the  ARCP  are  included  in  this  phase 
of  flight  planning.  Standard  EWO  planning  figures  were 
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Figure  3.  Flight  Planning  Phases 


used  for  the  time  and  distance  from  takeoff  to  level-off. 
These  figures  are  17  minutes  and  105  NM  for  the  bomber  and 
25  minutes  and  164  NM  for  the  tanker.  The  bomber's  level- 
off  gross  weight  is  476,000  pounds  while  the  tanker's 
level-off  gross  weight  varies  depending  on  its  takeoff 
gross  weight. 

After  level-off  both  the  bomber  and  tanker  follow 
a  maximum  range  cruise  (MRC)  profile  with  an  average  true 
airspeed  (TAS)  of  444  knots.  As  implied  by  its  name,  the 
MRC  profile  maximizes  aircraft  range  by  gradually  climbing 
the  aircraft  as  gross  weight  decreases  due  to  fuel  con¬ 
sumption.  The  ever  increasing  altitude  and  decreasing 
gross  weight  result  in  a  continually  decreasing  fuel  con¬ 
sumption  rate.  Since  they  are  constantly  changing,  the 
fuel  consumption  rate  and  altitude  are  recomputed  every  30 
minutes  in  accordance  with  performance  manual  procedures. 
The  net  result  of  these  computations  is  that  the  aircraft 
gross  weight  and  arrival  time  at  the  ARCP  can  be  computed. 
The  gross  weight  can  then  be  converted  to  fuel  load  by 
subtracting  the  aircraft  zero  fuel  weight.  These  standard 
EWO  weights  are  218,300  pounds  for  the  bomber  and  110,100 
pounds  for  the  tanker.  The  fuels  at  the  ARCP  are  important 
inputs  for  determining  onload  and  offload  capabilities. 

Holding  at  the  ARCP.  If  a  tanker  arrives  at  the 
ARCP  before  its  scheduled  bomber,  it  enters  a  holding 
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pattern  to  wait.  Holding  is  accomplished  at  maximum 
endurance  airspeed  to  minimize  fuel  consumption.  The  maxi¬ 
mum  endurance  airspeed  is  considerably  lower  than  cruise 
airspeeds  and  decreases  as  gross  weight  decreases.  Like 
the  cruise  portion  of  flight,  the  airspeed  and  fuel  con¬ 
sumption  rates  are  recomputed  every  30  minutes. 

Air  Refueling .  Air  refueling  commences  at  the 
ARCP  and  ends  at  the  end  air  refueling  (EAR)  point. 

Standard  planning  factors  include  an  altitude  of  30000 
feet,  TAS  of  400  knots,  and  a  fuel  transfer  rate  of  5000 
lbs/min  from  the  tanker  to  bomber. 

The  bomber's  capability  to  onload  fuel  is  predi¬ 
cated  on  completing  air  refueling  at  its  maximum  inflight 
gross  weight  of  488,000  pounds.  The  bomber's  average 
gross  weight  during  air  refueling  is  computed  as  follows: 

Average  Gross  Height  -  <89°°P-i--Start-tetuel^  "***  (l) 


The  average  gross  weight  is  then  used  to  compute  the  air 
refueling  fuel  consumption  rate.  This  fuel  consumption 
rate  in  lbs/min  is  subtracted  from  the  5000  lbs/min  transfer 
rate  to  obtain  a  net  transfer  rate.  The  time  required  f.or 
the  bomber  to  accomplish  air  refueling  is  then  determined 
by  the  following  equation. 


488000  -  Start  Gross  Weight 
Net  Transfer  Rate 


(2) 
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Multiplying  this  time  by  the  actual  transfer  rate  of 
5000  lbs/min  yields  the  bomber's  total  onload  capability. 
Its  net  onload  capability  is  488,000  minus  the  start  refuel 
ing  gross  weight  which  also  corresponds  to  the  total  onload 
capability  minus  the  fuel  consumed  while  obtaining  it. 

The  tanker’s  fuel  available  for  air  refueling  is 
determined  by  its  start  air  refueling  gross  weight  and  the 
weight  at  which  it  has  to  depart  for  its  PRB.  The  tanker's 
fuel  consumption  rate  during  refueling  is  also  based  on  its 
average  gross  weight  where 

Average  Cross  Weight  -  start  Q”38-  m  Weight  |3) 

Since  the  tanker  transfers  5000  lbs/min  of  fuel  to  the 
bomber,  it  has  a  net  transfer  rate  of  5000  lbs/min  plus 
its  fuel  bum  rate.  Dividing  the  fuel  available  for  air 
refueling  by  this  net  transfer  rate  gives  the  tanker's 
time  available  for  refueling. 

Start  Gross  Weight  -  ERR  Gross  Weight  ... 

Net  Transfer  Rate  m 

Its  total  offload  capability  is  the  product  of  this  time 
and  the  5000  lbs/min  transfer  rate.  The  remainder  of  the 
fuel  that  was  available  for  air  refueling  is  consumed  by 
the  tanker  itself. 
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Post  Air  Refueling .  After  completing  air  refuel¬ 
ing,  the  bomber  resumes  its  MRC  profile  to  either  the  next 
ARCP  or  to  the  entry  point,  whichever  is  applicable.  If 
it  is  the  former,  the  computations  just  discussed  under 
air  refueling  are  repeated.  If  it  is  the  latter  case, 
the  bomber's  entry  point  fuel  is  determined  by  subtracting 
its  zero  fuel  gross  weight  from  the  entry  point  gross 
weight. 

After  it  completes  refueling,  the  tanker  resumes 
the  MRC  profile  enroute  to  its  PRB .  The  tanker  has  to 
terminate  air  refueling  so  as  to  arrive  over  the  recovery 
base  with  its  required  fuel  reserve.  For  the  purposes  of 
this  study,  the  required  fuel  reserve  was  assumed  to  be 
5000  pounds  as  the  actual  figure  is  classified. 

Computing  the  tanker's  actual  EAR  gross  weight 
presents  a  slightly  different  problem  in  that  the  gross 
weight  over  the  PRB  is  known  and  the  EAR  wejight  is  unknown. 
This  problem  is  solved  by  flight  planning  in  reverse. 

Flight  planning  starts  over  the  recovery  base  and  proceeds 
backwards  in  30  minute  intervals  to  the  EAR  point. 

Air  Refueling  Location 

It  may  appear  from  Figure  2  and  the  discussion  to 
this  point  that  the  air  refueling  locations  are  fixed,  but 
this  is  not  the  case.  As  will  be  shown  below,  the  optimal 
location  of  the  ARCP  occurs  when  the  bomber's  total  onload 
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capability  equals  the  tanker's  total  offload  capability. 
These  capabilities,  in  turn,  are  functions  of  the  bomber 
and  tanker  fuel  loads  at  takeoff,  their  total  time  air¬ 
borne,  and  the  distance  to  the  tanker's  recovery  base. 

If  all  other  factors  are  held  constant,  the  bomber's  onload 
capability  increases  with  time  while  the  tanker's  offload 
capability  decreases.  These  are  the  natural  consequences 
of  fuel  consumption  over  time.  In  a  similar  manner,  the 
less  fuel  a  bomber  has  at  takeoff,  the  more  fuel  it  can 
onload  at  any  particular  time.  The  tanker,  on  the  other 
hand,  has  just  the  opposite  relationship.  The  less  fuel 
it  has  at  takeoff,  the  smaller  its  offload  capability  at 
any  particular  instant.  Finally,  as  the  distance  from  end 
air  refueling  to  the  tanker's  recovery  base  increases,  the 
offload  capability  decreases  because  of  the  additional 
fuel  required  to  reach  the  recovery  base. 

Bomber  onload  capability  and  tanker  offload  capa¬ 
bility  are  shown  as  functions  of  time  in  Figure  4.  For 
this  particular  combination  of  fuel  loads  and  recovery  base, 
offload  capability  equals  onload  capability  at  time  T*. 
Refueling  prior  to  this  time  will  result  in  the  tanker  land¬ 
ing  at  the  recovery  base  with  excess  fuel.  Refueling  after 
this  time  will  result  in  a  decreased  onload  for  the  bomber. 
The  time  T*  is  easily  equated  to  distance  and  thus  defines 
the  optimum  air  refueling  location  along  the  bomber's  route 
of  flight. 
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Bomber  Onload 


Figure  4.  Fuel  Onload-Off load  Capabilities 


Average  Tanker  Concept 


Since  all  of  the  bombers  are  assumed  to  takeoff 
at  the  same  gross  weight,  their  onload  capability  is 
strictly  a  function  of  time.  This  is  not  the  case  for  the 
tankers.  Their  takeoff  fuel  loads  can  vary,  and  this 
affects  their  offload  capability.  Their  offload  capabil¬ 
ity  is  also  affected  by  distance  to  the  PRB.  As  a  result, 
it  became  necessary  to  define  an  average  tanker  for  each 
of  the  three  possible  air  refuelings.  As  will  be  seen  in 
Chapter  III,  these  average  tankers  are'  used  to  determine 
the  number  of  tankers  required  by  each  bomber  and  to  estab¬ 
lish  initial  solutions  for  the  mating  problem. 

The  average  tankers  were  defined  by  first  computing 
the  bomber's  onload  capability  in  30-minute  increments 
starting  at  30  minutes  after  level-off.  The  tanker's  off¬ 
load  capability  was  computed  over  the  same  increments  based 
on  a  maximum  takeoff  fuel  load  and  PRBs  that  were  one,  two, 
and  three  hours  from  the  EAR  point.  These  times  to  the 
PRB  were  selected  as  representative  since  the  actual  time 
might  vary  anywhere  from  just  a  few  minutes  to  three  hours 
or  more.  The  resulting  onload  and  offload  capabilities 
are  plotted  graphically  in  Figure  5.  The  optimum  refueling 
points  range  from  approximately  3  hours  and  55  minutes  to 
4  hours  and  20  minutes  after  takeoff  depending  on  the 
PRB.  As  is  to  be  expected,  a  tanker  with  one  hour  PRB  has 
the  highest  offload  capability  and  a  tanker  with  the  three 


hour  PRB  has  the  lowest  offload  capability.  A  complete 
listing  of  the  parameters  for  the  first  refueling  is  shown 
in  Table  1. 


TABLE  Z 

FIRST  REFUELING  PARAMETERS 


Tims  to 
PRB 

Time  from 
Takeoff 
to  ARCP 

Distance  from 
Takeoff  to 
ARCP 

Qnload/Off- 
load  (1000s 
of  lbs) 

Time  on 
Track 

Track 

Length 

1  hr 

4.33  hrs 

1887  NM 

100.5 

20.1  min 

134  m 

2  hr 

4.13  hrs 

1800  NM 

96.0 

19.0  min 

130  NM 

3  hr 

3.90  hrs 

1697  NM 

92.0 

18.4  min 

123  NM 

A  bomber's  onload  capability  at  any  given  point  for 
the  second  refueling  depends  on  its  first  refueling.  If 
its  first  refueling  was  with  a  tanker  that  had  a  distant 
PRB,  its  onload  capability  will  be  more  than  if  it  refueled 
with  a  tanker  that  had  a  close-in  PRB.  This  is  reflected 
by  three  onload  curves  for  the  second  refueling  of  Figure  5. 
These  onload  curves  correspond  to  the  three  possibilities 
for  the  first  refueling.  They  are  identical  to  the  original 
onload  curve  except  they  have  been  shifted  to  the  right  by 
the  amount  it  took  to  complete  the  first  refueling,  and  by 
the  time  it  takes  the  bomber  to  reduce  its  gross  weight  to 
the  same  weight  it  had  30  minutes  after  initial  level-off. 
These  three  onload  curves  then  combine  with  the  three 
possible  tankers  to  produce  nine  onload  curves  for  the 


third  air  refueling;  however,  only  the  onload  curves  cor¬ 
responding  to  the  one,  two,  and  three  hour  PRBs  of  the 
second  refueling  are  shown.  Once  again  these  curves  reflect 
the  original  curve  shifted  by  the  appropriate  times. 

Out  of  all  of  these  possible  air  refueling  combina¬ 
tions,  the  tanker  with  the  two  hour  recovery  base  was 
selected  as  the  average.  These  combinations  and  the  optimum 
refueling  times  are  denoted  by  the  dashed  lines  of  Figure  5. 
The  characteristics  of  these  average  tankers  are  listed 
in  Table  II. 


TABLE  II 

AVERAGE  TANKER  CHARACTERISTICS 


Refuel¬ 

ing 

* 

Time  from 
Takeoff 
to  ARCP 

Distance  from 
Takeoff  to 
ARCP 

Onload/Off¬ 
load  (1000s 
of  lbs) 

Time  on 
Track 

Track 

Length 

1 

4.13  hrs 

1880  rat 

96.0 

19.0  min 

130  NM 

2 

7.23  hrs 

3160  NM 

63.0 

12.8  min 

85  NM 

3 

9.27  hrs 

4060  m 

44.0 

9.0  min 

60  NM 

Variations  in  the  tanker  offload  capabilities 
required  the  development  of  the  average  tanker.  The  average 
tanker,  in  turn,  requires  procedures  for  dealing  with  devi¬ 
ations  from  the  average.  The  result  of  these  deviations  is 
that  the  offload  capability  will  no  longer  equal  the  onload 
capability.  This  can  be  compensated  for  by  shifting  the 
ARCP.  If  the  offload  capability  is  greater  than  the  onload 
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capability,  the  distance  to  the  ARCP  has  to  be  increased. 

Such  a  shift  decreases  the  offload  capability  and  increases 
the  onload  capability  because  more  fuel  is  required  to 
reach  the  ARCP.  When  the  offload  capability  is  less  than 
the  onload  capability,  the  distance  to  the  ARCP  has  to  be 
decreased.  This  shift  increases  the  offload  capability 
and  decreases  the  onload  capability  because  less  fuel  is 
consumed  to  reach  the  ARCP.  If  the  bomber  has  more  than 
one  refueling,  shifting  the  first  ARCP  will  require  a  shift 
in  the  second  ARCP  which  will  then  require  a  shift  in  the 
third  ARCP,  if  applicable.  Shifting  each  of  these  ARCPs 
the  appropriate  distances  will  equalize  offload  and  onload 
capabilities. 

There  is  one  final  and  important  point  that  arises 
from  the  average  tanker  concept.  Mating  one  bomber  with 
one  tanker  that  could  go  to  three  PRBs  generated  three 
refueling  combinations  for  the  first  refueling,  nine  for 
the  second,  and  twenty- seven  for  the  third.  Since  the  time 
for  any  tanker  to  reach  any  PRB  will  rarely  (if  ever )  be 
identical,  the  end  result  is  that  the  total  number  of  refuel¬ 
ing  combinations  for  any  given  refueling  is  the  product  of 
the  number  of  tankers  and  PRBs  available.  For  example, 
assume  that  there  are  ten  tankers  and  recovery  bases  avail¬ 
able.  This  would  result  in  up  to  100  possible  refueling 
combinations  for  a  bomber  on  the  first  refueling.  This 
would  generate  100  onload  curves  for  the  second  refueling. 
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each  of  which  can  again  combine  with  100  possible  tanker 
and  PRB  combinations.  This  would  result  in  10000  possible 
alternatives,  each  of  which  can  again  combine  with  100 
possible  tanker  and  PRB  combinations.  Thus,  ten  tankers 
and  PRBs  can  generate  up  to  one  million  refueling  combina¬ 
tions  for  only  one  bomber!  Expanding  the  problem  to  several 
hundred  bombers,  tankers,  and  PRBs  would  generate  an 
incredibly  large  number  of  refueling  combinations.  This 
eliminates  enumeration  as  an  effective  method  of  optimizing 
the  mating  process. 

Computerized  Flight  Planning 

All  of  the  flight  planning  to  this  point  has  been 
accomplished  by  manual  look-up  in  the  appropriate  aircraft 
performance  manuals.  Such  an  approach  is  obviously  not 
amenable  to  developing  computerized  algorithms  for  solving 
the  mating  problem;  therefore,  an  attempt  was  made  to 
obtain  the  computerized  performance  polynomials  used  by 
SAC  and  MARP.  Obtaining  these  programs  turned  out  to  be 
difficult,  and  there  was  no  guarantee  that  they  could  be 
adapted  to  the  CDC  computers  if  they  were  obtained.  For 
these  reasons,  linear  regression  techniques  were  used  to 
develop  equations  for  the  appropriate  performance  parameters. 
These  parameters  included: 

1.  B-52  maximum  range  cruise  fuel  consumption 

2.  B-52  fuel  consumption  during  air  refueling 


at  30000  feet 


3.  KC-135  maximum  range  cruise  fuel  consumption 

4.  KC-135  fuel  consumption  during  maximum 
endurance  holding  at  30000  feet 

5.  KC-135  fuel  consumption  during  air  refueling 
at  30000  feet 

All  of  these  parameters  are  functions  of  gross 
weight  and  altitude.  For  the  constant  altitude  conditions 
and  small  gross  weight  ranges  of  parameters  2,  4,  and  5, 
fuel  consumption  varies  almost  linearly  with  the  gross 
weight.  This  is  not  the  case  for  parameters  1  and  3.  The 
altitude  is  not  constant,  and  the  gross  weight  varies  over 
a  wider  range;  however,  these  variations  can  be  accounted 
for  by  dividing  the  gross  weight  range  into  two  smaller 
ranges.  When  this  is  done,  fuel  consumption  again  varies 
almost  linearly  with  gross  weight  over  each  of  these 
smaller  ranges. 

Linear  regressions  were  run  for  fuel  flow  in  pounds/ 
minute  versus  gross  weight  in  thousands  of  pounds  for 
parameters  1  through  5  using  the  Statistical  Package  for 
the  Social  Sciences  (SPSS)  routines.  The  input  data,  the 
resulting  equations,  and  a  summary  of  the  SPSS  printouts 
are  included  in  Appendix  A.  The  correlation  coefficient 
of  each  of  these  equations  exceeded  .98  which  indicates 
the  high  degree  of  linearity  between  gross  weight  and  fuel 
flow.  These  high  correlation  coefficients  combined  with 
the  fact  that  these  equations  are  used  for  all  three 
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approaches  are  deemed  as  ample  justification  for  using 
the  linear  regression  results  in  lieu  of  the  performance 
polynomials. 


The  second  problem  encountered  in  the  flight  plan¬ 
ning  process  was  computation  of  the  great  circle  distances. 
This  is  accomplished  by  using  two  subroutines  adapted  from 
Reference  1.  Subroutine  Circle  computes  the  great  circle 
distance  between  two  points  when  given  the  coordinates  of 
the  points.  Subroutine  Latlon  yields  the  coordinates  of 
a  second  point  given  the  coordinates  of  the  first  point 
and  the  great  circle  distance  and  course.  These  two  sub¬ 
routines  are  included  in  Appendix  C. 

Summary 

The  prototype  problem,  flight  planning  process, 
and  the  concepts  of  cost  and  average  tanker  were  introduced 
and  developed  in  this  chapter  because  these  items  are 
pertinent  to  the  discussion  of  theory  and  methodology 
which  follows  in  Chapter  III.  Their  treatment  at  this 
point  provides  the  necessary  background  information  for 
this  discussion. 
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III.  Theory  and  Methodology 

Two  new  methods  are  explored  for  solving  the 
bomber  and  tanker  mating  problems.  They  are  a  network 
method  and  a  "greedy"  method.  The  first  method  uses  net¬ 
work  theory  in  an  attempt  to  obtain  an  optimal  or  near 
optimal  assignment  of  bombers,  tankers,  and  recovery  bases. 
The  second  method  uses  a  marginal  cost  improvement  algorithm 
to  make  these  assignments.  Neither  method  offers  a  guar¬ 
antee  of  optimality;  however,  the  second  is  easy  to  imple¬ 
ment  and  similar  "greedy”  algorithms  have  been  employed  in 
a  wide  variety  of  applications  with  varying  degrees  of  suc¬ 
cess.  Thus,  it  serves  as  a  basis  of  comparison  for  both 
the  network  method  and  the  current  Logicon  method. 

The  remainder  of  this  chapter  is  devoted  to  the 
theory  and  methodology  behind  these  two  methods  as  well 
as  the  Logicon  method.  The  network  method  is  discussed 
first.  It  is  followed  by  a  discussion  of  the  Logicon 
method.  The  Logicon  method  is  included  in  this  discussion 
because  it  is  necessary  to  understand  how  it  works  in  order 
to  be  able  to  better  compare  it  with  the  network  method. 
Finally,  the  theory  and  methodology  of  the  "greedy"  method 
are  discussed. 
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Network  Method 


The  underlying  concept  of  the  network  method  is 
the  network.  Network  models  have  been  used  to  solve  a 
variety  of  very  complex  problems  that  include,  but  are  not 
limited  to  transportation  of  goods,  design  of  communication 
and  pipeline  systems,  assignment  of  men  to  jobs,  bid  evalua¬ 
tion,  and  production  planning  (Ref  2:1). 

According  to  the  terminology  of  the  theory  of 
graphs,  a  graph  consists  of  a  set  of  junction  points  called 
nodes,  with  certain  pairs  of  the  nodes  being  joined  by 
lines  called  arcs  (Ref  10:234)  .  Figure  6  is  an  example  of 
a  graph  where  the  circles  are  the  nodes.  They  are  desig¬ 
nated  as  1,  2,  3,  and  4.  These  nodes  are  connected  by  the 
arcs  (1,2),  (1,3),  (2,3),  (3,2),  (2,4),  and  (3,4).  As  can 
be  seen  in  this  example,  all  nodes  do  not  have  to  be  con¬ 
nected,  e.g.,  nodes  1  and  4  are  not  connected. 


A  network  is  a  graph  with  flow  in  its  arcs,  and  is 
said  to  be  directed  if  its  arcs  are  oriented  in  a  specified 
direction.  If  the  arcs  in  Figure  6  had  flow  in  them,  this 
figure  would  be  an  example  of  a  directed  network  with  the 
flow  indicated  by  the  directional  arrows  on  the  arcs. 

The  approach  of  the  network  method  is  to  formulate 
the  bomber  and  tanker  mating  problem  as  a  directed  network 
problem.  The  nodes  of  the  network  are  the  tanker  bases, 
ARCPs,  EAR  points,  and  PRBs.  The  tankers  "flow"  through 
arcs  from  their  bases  to  the  ARCPs,  EAR  points,  and  PRBs, 
respectively.  The  objective  is  to  flow  the  tankers  through 
this  network  at  the  minimum  cost  which  should  maximize 
bomber  entry  point  fuels.  As  indicated  in  the  previous 
chapter  these  costs  represent  the  tanker  fuel  used  in 
traversing  the  arcs. 

More  specifically,  the  network  method  is  structured 
similar  to  the  capacitated  transshipment  model  (also  known 
as  the  minimum  cost  flow  problem)  which  determines  in  what 
quantities  or  at  what  rates  a  good  should  flow  through  the 
arcs  of  a  network  so  as  to  minimize  total  shipment  costs 
(Ref  6:3).  The  arcs  of  the  network  consist  of  ordered 
pairs  of  nodes  (tail  to  head)  and  are  indexed  by  k.  Each 
arc  has  a  shipping  cost  per  unit  of  flow,  C^,  a  minimum 
allowable  flow  (lower  bound),  L^,  and  a  maximum  allowable 
flow  (upper  bound),  U^.  The  nodes  of  the  network  are  either 
supply  nodes  where  units  enter  the  network,  demand  nodes 


where  units  depart  the  network,  or  transshipment  nodes 
where  the  units  just  pass  through.  The  capacitated  trans¬ 
shipment  model  minimizes  the  total  costs  with  flows  that 
satisfy  the  associate  upper  and  lower  bounds  and  preserve 
the  conservation  of  flow  at  each  node.  Mathematically, 
this  can  be  expressed  as 

Minimize:  k  e  A  CkXk 

Subject  to:  k  *  k  k  e  A  *k  *  bi  for  1  £  “ 

with  tail  i  with  head  i 

Lk  1  \  i  uk  for  *  e  A 

where  b^  *  Supply  if  i  is  a  supply  node 
-Demand  if  i  is  a  demand  node 
0  otherwise 

and  A  is  the  set  of  all  arcs 

N  is  the  set  of  all  nodes 

This  problem  can  be  solved  by  using  linear  program¬ 
ming  techniques  or  by  using  one  of  several  special  purpose 
network-flow  computer  programs.  These  latter  programs  can 
solve  these  problems  up  to  200  times  faster  than  most 
typical  linear  programming  codes  by  taking  advantage  of 
the  special  network  structure.  One  such  program  is  GNET, 
and  it  has  been  incorporated  into  the  network  method  as  a 
subroutine.  It  uses  a  primal-simplex  method  to  solve  the 
capacitated  transshipment  problem.  This  approach  is  con¬ 
sidered  to  be  much  more  efficient  than  most  of  the  other 
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programs  which  generally  use  an  out-of -kilter  approach 
(Ref  2:3) .  A  full  description  of  GNET  and  its  capabili¬ 
ties  is  contained  in  Appendix  D. 

Initial  Network  Method.  The  methodology  of  the 
initial  network  method  is  outlined  in  Figure  7.  It  is 
referred  to  as  the  initial  network  method  because  it  later 
turned  out  to  be  infeasible,  and  had  to  be  altered  slightly 
Each  of  the  blocks  or  steps  of  this  method  are  discussed 
in  turn  using  the  prototype  problem  described  in  Chapter  II 
for  illustration  purposes. 

The  first  step  consists  of  determining  how  many 
tankers  should  be  assigned  to  each  bomber.  This  is  done 
by  flying  the  bomber  to  its  entry  point  unrefueled  and 
noting  how  much  fuel  it  arrives  with.  This  figure  is  com¬ 
pared  with  the  required  entry  point  fuel  to  determine  the 
additional  fuel  required.  The  additional  fuel  required  is 
then  equated  to  the  number  of  tankers  needed.  Since  the 
offload  capability  of  a  tanker  depends  on  its  takeoff 
gross  weight  and  recovery  base,  the  average  tanker  of 
Chapter  II  is  used  to  make  this  determination.  For  con¬ 
venience,  the  average  tanker  offload  capabilities  are 
repeated  in  Table  III;  however,  there  is  one  final  adjust¬ 
ment  that  has  to  be  made  to  these  figures.  As  noted  in 
Chapter  II,  the  bomber  consumes  part  of  its  onload  in  the 
process  of  obtaining  it.  In  addition,  the  bomber  is 


Determine  Number  of  Tankers 
to  be  Assigned  to  Each  Bomber 


Figure  7.  initial  Methodology  for  the  Network  Method 


TABLE  III 


OFFLOAD  CAPABILITIES  IN  THOUSANDS  OF  POUNDS 


Average  Offload  Effective  Offload  Cumulative  Effective 
Refueling  Capability  Capability  Offload  Capability 


1 

96.0 

83.0 

83.0 

2 

63.0 

52.0 

135.0 

3 

44.0 

29.0 

164.0 

heavier  after  refueling,  and  thus  burns  more  fuel  than  it 
would  if  unrefueled.  The  net  result  of  these  two  factors 
is  that  the  bomber's  entry  point  fuel  with  refueling  is 
less  than  the  sum  of  the  tanker's  offload  capability  and 
the  bomber's  entry  point  fuel  without  the  refueling ( s) . 

This  is  reflected  in  the  effective  offload  figures  of 
Table  III.  These  figures  were  found  by  computing  the 
average  net  gain  in  entry  point  fuel  for  numerous  refueling 
situations.  The  last  column  of  Table  III  consists  of  the 
cumulative  effective  offload  for  one,  two,  and  three  refuel¬ 
ings.  These  are  the  numbers  that  are  used  to  determine 
how  many  tankers  are  actually  required  to  meet  the  bomber's 
entry  point  fuel  requirement. 

Table  IV  shows  the  entry  point  fuels  without  refuel¬ 
ing,  the  required  entry  point  fuels,  and  the  additional 
fuel  required  for  each  of  the  four  bombers  of  the  prototype 
problem.  Also  shown  are  the  tankers  that  would  be  required 
to  supply  the  additional  fuel  requirements.  If  all  of 
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TABLE  IV 


PROTOTYPE  PROBLEM  FUEL  REQUIREMENTS  IN 
THOUSANDS  OF  POUNDS 


Bcmber 

Entry  Point 
Fuel 

unrefueled 

Entry  Point 
Fuel 
Required 

Difference 

Tankers 

Required 

Surplus 

1 

123.7 

205.0 

81.3 

1 

1.7 

2 

90.5 

250.0 

159.5 

3 

4.5 

3 

118.1 

220.0 

101.9 

2 

33.1 

4 

98.8 

250.0 

151.2 

3 

12.8 

these  tankers  were  available,  each  bomber  would  receive 
fuel  in  excess  of  its  requirements.  The  excess  can  be  esti- 
mated  by  subtracting  the  required  entry  point  fuel  from 
the  sum  of  the  cumulative  effective  offload  and  the  unre¬ 
fueled  entry  point  fuel.  These  figures  are  shown  in  the 
last  column  of  Table  IV. 

As  can  be  seen  from  Table  IV,  nine  tankers  are 
required  to  meet  all  of  the  bombers '  entry  point  fuel 
requirements,  but  there  are  only  seven  tankers  available 
in  the  prototype  problem.  Two  refuelings  have  to  be 
deleted.  This  is  accomplished  on  the  basis  of  which  bomber 
has  the  greatest  fuel  surplus.  This  procedure  deletes 
refuelings  for  those  bombers  that  can  best  afford  it.  This 
process  is  continued  until  the  tankers  required  equal  the 
tankers  available.  For  the  prototype  problem,  bombers  3 


and  4  have  the  largest  surpluses.  Each  of  them  lose  one 
refueling  respectively. 

Once  the  number  of  tankers  required  for  each 
bomber  has  been  determined,  the  network  can  be  defined. 

The  network  formulation  of  the  prototype  problem  is  shown 
in  Figure  8 .  The  source  node  starts  the  flow  of  tankers 
to  the  tanker  base  nodes.  These  nodes  are  then  connected 
to  each  refueling  track,  and  from  there  to  each  post¬ 
refueling  base.  The  post-refueling  bases  are  in  turn  con¬ 
nected  to  the  secondary  sink  node,  and  the  secondary  sink 
node  is  connected  to  the  primary  sink  node.  Only  one  sink 
node  is  required  for  a  network  problem  in  general,  but  two 
were  required  in  this  formulation  because  of  Subroutine 
GNET.  It  requires  two  sinks  because  the  total  supply  from 
the  source  node  must  equal  the  total  demand  at  the  sink 
node.  Since  the  sum  of  all  the  PRB  capacities  exceeds  the 
total  supply  of  tankers  that  emanated  from  the  source 
node,  the  primary  sink  node  is  required  to  equate  the  total 
number  of  tankers  initially  available  to  the  number  that 
flow  into  the  primary  sink.  For  the  sake  of  clarity,  only 
the  arcs  connecting  tanker  base  1  with  each  refueling  track 
and  refueling  track  1  with  each  post-refueling  base  are 
shown.  In  reality,  every  tanker  base  is  connected  to  each 
refueling  track,  and  every  refueling  track  is  connected  to 
each  post-refueling  base. 
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The  next  step  of  the  network  method  involves  set¬ 
ting  the  lower  and  upper  bounds  for  the  flows  through  the 
various  arcs  in  the  network.  Each  arc  from  the  source 
node  to  a  tanker  base  node  has  a  lower  bound  of  zero  and  an 
upper  bound  equal  to  the  number  of  tankers  assigned  to  that 
base.  The  arcs  from  the  tanker  base  nodes  to  the  ARCP 
nodes  have  a  lower  bound  of  zero  and  an  upper  bound  of  one. 
The  lower  and  upper  bounds  of  the  air  refueling  arcs  are 
both  one.  This  insures  that  an  air  refueling  takes  place 
on  each  track.  The  arcs  from  the  EAR  points  to  the  PRBs 
have  a  lower  bound  of  zero  and  an  upper  bound  of  one.  The 
lower  bounds  of  the  arcs  from  the  PRB  to  the  secondary  sink 
are  zero  and  the  upper  bounds  are  equal  to  each  recovery 
base's  capacity.  Finally,  the  arc  into  the  primary  sink 
has  a  lower  bound  of  zero  and  an  upper  bound  equal  to  the 
total  number  of  tankers  available.  As  mentioned  previously, 
this  last  arc  insures  that  supply  equals  demand.  The  respec¬ 
tive  upper  and  lower  bounds  for  each  arc  are  enclosed  in 
parentheses  in  Figure  8 . 

The  fourth  step  of  the  network  method  computes  the 
costs  of  flowing  the  tankers  through  the  arcs.  In  this 
particular  application,  these  costs  represent  the  fuel  con¬ 
sumed  while  traversing  the  arcs.  The  costs  of  flowing  the 
tankers  from  the  source  to  the  tanker  bases  and  from  the 
recovery  bases  to  the  sinks  are  zero  because  these  arcs 
are  only  required  to  establish  the  flow.  They  do  not 


affect  the  total  cost  function.  The  remainder  of  the  costs 
are  computed  using  the  flight  planning  process  detailed  in 
Chapter  II.  These  costs  include  the  differences  in  tanker 
takeoff  gross  weights,  the  fuel  consumed  from  takeoff  to 
start  refueling  (including  holding  if  applicable),  and 
the  fuel  consumed  from  the  EAR  point  to  the  PRB.  The  first 
two  costs  are  aggregated  into  one  figure  by  subtracting  the 
fuel  consumed  to  reach  the  ARCP  from  the  takeoff  gross 
weight.  This  figure  takes  into  account  the  fact  that, 
although  a  lighter  tanker  consumes  less  fuel  to  reach  the 
ARCP,  it  will  still  have  less  offload  capability  than  a 
heavier  tanker.  Unlike  the  other  costs,  it  is  obvious  that 
this  figure  should  be  maximized  in  order  to  maximize  the 
offload.  GNET,  on  the  other  hand,  attempts  to  minimize 
costs.  This  discrepancy  is  overcome  by  defining  this 
refueling  cost  as  a  negative  cost.  Thus,  minimizing  the 
negative  cost  is  equivalent  to  maximizing  the  fuel  avail¬ 
able  at  the  ARCP.  One  final  consideration  concerns  infeas¬ 
ible  bomber  and  tanker  matings.  The  refueling  costs  for 
these  assignments  are  set  at  a  very  large  positive  number. 
This  prevents  these  tankers  from  being  considered  in  the 
final  solution.  These  arcs  could  have  also  been  eliminated 
from  the  network,  but  were  retained  for  ease  of  computeriza¬ 
tion. 

The  final  step  of  the  network,  method  consists  of 
solving  for  the  optimal  flow  of  tankers  through  the  network 
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Unfortunately,  this  last  step  proved  to  be  impossible  for 
this  particular  formulation  of  the  problem.  To  obtain  the 
optimal  mating  of  bombers  and  tankers,  it  is  necessary  to 
assign  the  optimal  tanker  to  the  optimal  refueling  location 
and  the  optimal  post-refueling  base  for  each  possible 
refueling;  however,  the  optimal  refueling  location  is  a 
function  of  the  tanker  and  recovery  base  assignments  which 
are  in  turn  functions  of  the  refueling  location.  This  type 
of  problem  is  referred  to  as  a  three-dimensional  assignment 
problem,  and  belongs  to  a  class  of  problems  known  as 
NP-complete  problems.  There  is  no  known  polynomially 
bounded  algorithm  that  is  able  to  solve  problems  in  this 
class  (Ref  14:8-9).  This  obviously  required  a  reformula¬ 
tion  of  the  problem. 

The  Revised  Network  Method.  Since  there  is  no 
efficient  procedure  for  solving  the  mating  problem  described 
above,  an  alternative  or  heuristic  approach  is  required. 

The  revised  network  method  is  such  an  approach.  This 
approach  utilizes  the  second  aspect  of  the  average  tanker 
concept  of  Chapter  II  to  eliminate  one  dimension  of  the 
three-dimensional  assignment  problem.  The  solution  that  is 
obtained  is  then  iterated  in  an  attempt  to  further  improve 
the  solution. 

The  revised  network  method  is  outlined  in  Figure  9. 
The  only  change  in  the  first  four  steps  from  those  of  the 
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Determine  the  Number  of  Tankers 
to  be  Assigned  to  Each  Bomber 


Figure  9.  Revised  Methodology  for  the  Network  Method 


initial  method  is  the  manner  in  which  the  network  is  defined. 
The  basic  structure  as  shown  in  Figure  8  is  unchanged;  the 
only  difference  is  that  the  ARCP  and  EAR  point  locations 
are  no  longer  dependent  on  the  tanker  and  PRB.  They  are 
fixed  as  explained  in  the  next  paragraph. 

One  dimension  of  the  three-dimensional  assignment 
problem  of  step  5  is  eliminated  by  assuming  that  average 
tankers  are  assigned  to  each  air  refueling.  This  assump¬ 
tion  fixes  the  locations  of  the  refueling  tracks.  The  dis¬ 
tances  from  the  bomber's  departure  base  to  the  ARCPs  were 
calculated  for  average  tankers  in  Chapter  II.  These  dis¬ 
tances,  the  corresponding  track  lengths,  and  end  air  refuel¬ 
ing  distances  are  repeated  in  Table  V.  Since  the  air 
refueling  location  is  fixed  and  no  longer  dependent  on  the 
tanker  and  PRB  assignments,  the  GNET  subroutine  is  able  to 
flow  the  tankers  through  the  network.  The  net  result  is 
that  each  tanker  is  assigned  to  the  refueling  tracks  and 
PRBs  so  as  to  minimize  the  total  cost.  This  provides  an 
initial  solution  to  the  mating  problem. 

In  reality  however,  the  actual  tanker  assigned  to 
a  refueling  track  is  rarely  an  average  tanker.  This  means 
the  refueling  is  not  optimal  because  the  offload  and  onload 
capabilities  are  not  equal.  The  sixth  step  of  the  revised 
network  method  optimizes  these  refuelings.  It  does  this 
by  adjusting  the  air  refueling  locations  (as  described  in 
Chapter  II)  until  the  offload  and  onload  capabilities  are 
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TABLE  V 


REFUELING  LOCATIONS 

FOR  AVERAGE 

TANKERS 

Refueling 

Distance  from 
Takeoff  to  ARCP 

Track 

Length 

Distance  from 
Takeoff  to  EAR 

1 

1800  m 

130  1*4 

1930  NM 

2 

3160  NM 

85  NM 

3245  NM 

3 

4060  NM 

60  NM 

4120  NM 

within  400  pounds  of  each  other.  This  400-pound  tolerance 
is  a  compromise  between  the  desired  accuracy  and  the  com¬ 
puter  time  required  to  achieve  it.  The  latter  becomes  a 
factor  in  large  problems  with  many  multiple  refuelings 
because  adjusting  the  first  refueling  location  also 
requires  adjusting  the  second  and  third  refueling  locations 
as  applicable. 

These  adjustments  to  the  air  refueling  locations 
also  affect  the  arc  costs  which  determined  the  assignments 
to  begin  with.  The  last  step  of  the  revised  network  method 
iterates  the  assignment  process  in  an  attempt  to  improve 
the  solution.  Each  successive  iteration  uses  the  refueling 
locations  from  the  previous  iteration  as  the  new  initial 
solution.  The  number  of  iterations  desired  is  determined 
by  the  user.  For  the  purposes  of  this  study,  this  number 
was  initially  set  at  ten. 

The  revised  network  method  was  settled  upon  as  the 
network  approach  for  solving  the  mating  problem.  It  will 


thus  be  referred  to  as  the  network  method  for  the  remainder 
of  this  report. 

Current  Method 

Since  the  objective  of  this  research  effort  is  to 
develop  an  improved  method  for  solving  the  bomber  and 
tanker  mating  problem,  it  is  desirable  to  compare  any  pro¬ 
posed  methods  to  the  one  currently  in  use.  This  method  is 
Logicon  Corporation’s  Mating  and  Ranging  Program  (MARP) . 

It  is  an  extremely  large  and  complex  program  that  performs 
many  other  functions  in  addition  to  solving  the  mating 
problem.  It  is  also  written  in  an  advanced-language  that 
is  incompatible  with  the  AFIT  computers.  As  a  result,  it 
was  not  possible  to  use  MARP  itself  in  this  study.  Instead 
a  separate  program  was  developed  that  emulates  the  method¬ 
ology  used  by  MARP  in  the  assignment  process.  This  program 
referred  to  as  the  pseudo-Logicon  method,  is  then  used  as 
the  basis  of  comparison.  The  methodologies  of  MARP  and  the 
pseudo-Logi con  method  which  was  developed  for  this  study 
are  discussed  below. 

MARP.  The  methodology  of  MARP  is  similar  to  that 
of  the  network  method  in  that  it  also  uses  network  theory 
to  obtain  an  optimal  or  near  optimal  solution  to  the^-bomber 
and  tanker  mating  problem.  It  also  defines  costs  in  a 
manner  similar  to  the  network  method,  and  then  flows  the 
tankers  through  a  network  to  minimize  these  costs.  The 


network  has  a  slightly  different  structure  because  a  dif¬ 
ferent  network  solving  algorithm  is  used.  This  algorithm 
is  known  as  PNET  and  also  uses  a  primal -simplex  method  to 
solve  the  capacitated  transshipment  problem  (Ref  12) . 

Up  to  this  point,  there  are  very  few  differences 
between  MARP  and  the  network  method.  There  is,  however, 
one  major  difference.  This  concerns  their  handling  of  the 
post-refueling  bases.  MARP  assigns  each  tanker  to  the  best 
(closest)  post-refueling  base  without  regard  to  the  recovery 
base's  capacity.  After  all  assignments  have  been  made,  it 
then  checks  to  see  if  any  PRB  capacities  have  been  exceeded. 
If  they  have  been,  it  reassigns  the  excess  tankers  to 
other  unsaturated  PRBs.  These  tankers  are  reassigned  on 
the  basis  of  their  bombers'  entry  point  fuel  states.  First, 
they  are  ranked  in  the  order  of  weakest  to  strongest  entry 
point  fuel  state  where  the  weakest  bomber  is  the  one  that 
is  the  furthest  below  its  desired  entry  point  fuel.  Then 
the  tanker  associated  with  the  weakest  state  is  reassigned 
first.  It  is  sent  to  the  next  best  (closest)  PRB  relative 
to  its  EAR  point.  This  process  is  continued  until  all 
tankers  are  reassigned.  Assigning  them  from  the  weakest 
to  the  strongest  insures  that  the  tankers  most  able  to 
afford  it  are  reassigned  to  the  farthest  PRBs  (Ref  15) . 

Since  changing  the  PRBs  affects  the  tankers'  off¬ 
load  capabilities,  MARP  then  readjusts  the  refueling  loca¬ 
tions  so  as  to  equalize  offload  and  onload  capabilities. 


The  final  step  of  the  MARP  method  perturbates  the 
solution  to  see  if  it  can  be  improved.  This  step  consists 
of  aribitrarily  changing  a  limited  number  of  tanker  assign¬ 
ments  and  seeing  if  any  improvements  are  made  in  the  entry 
point  fuels.  If  an  improvement  is  attained,  these  new 
tanker  assignments  become  the  final  solution.  Otherwise, 
the  original  solution  stands. 

P seudo -Log ic on  Method .  This  method,  as  developed 
by  the  authors,  duplicates  MARP  through  a  three-step  pro¬ 
cess.  The  first  step  uses  the  network  method  to  obtain  the 
initial  bomber  and  tanker  matings  without  regard  to  PRB 
capacities.  This  is  accomplished  by  making  each  PRB 
capacity  equal  to  or  greater  than  the  total  number  of 
tankers  available. 

The  second  step,  like  MARP,  checks  each  PRB  to  see 
if  its  actual  capacity  has  been  exceeded.  If  so,  it 
reassigns  the  excess  tankers  using  the  same  logic  as  MARP 
and  readjusts  the  refueling  locations. 

The  final  step  then  repeats  or  iterates  the  entire 
process  as  in  the  network  method.  That  is,  the  refueling 
locations  of  the  initial  solution  become  the  new  fixed 
locations  for  the  next  iteration  of  the  network  method. 

This  step,  if  anything,  should  be  superior  to  the  limited 
perturbations  of  the  MARP  method. 
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It  is  not  claimed  that  the  pseudo-Logicon  method 
is  identical  to  the  MARP  method,  but  is  believed  to  be 
close  enough  to  serve  as  a  basis  of  comparison  with  other 
methods.  In  fact,  the  iterative  process  may  be  an  improve¬ 
ment  over  MARP.  If  this  is  the  case,  it  biases  the  com¬ 
parisons  in  favor  of  MARP. 

The  "Greedy"  Method 

The  "greedy"  method  is  similar  to  the  Vogel  Approxi¬ 
mation  Method  which  has  enjoyed  widespread  use  as  a  method 
of  finding  an  initial  feasible  solution  to  a  transportation 
problem  (Ref  10:134).  The  "greedy"  method  developed  for 
this  study  employs  a  fairly  simple  algorithm  and  is  best 
illustrated  with  an  example.  Such  an  example  is  shown  in 
Figure  10.  This  example  is  a  typical  assignment  problem 
where  the  objective  is  to  assign  the  machines  to  the  jobs 
at  minimum  cost. 

The  first  step  is  to  compute  the  costs,  if  neces¬ 
sary.  For  this  example,  the  costs  are  given.  The  next 
step  is  to  find  the  difference  between  the  smallest  and 
next  smallest  cost  in  each  row.  These  differences  are 
shown  in  the  difference  column  of  Figure  10.  The  first 
machine  to  be  assigned  is  the  one  with  the  largest  differ¬ 
ence.  This  corresponds  to  machine  4  of  Figure  10a. 

Machine  4  is  then  assigned  to  the  job  that  results  in  the 
lowest  cost.  This  is  job  3  which  is  circled.  Machine  4 
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Machine 


Job 


Difference 


and  job  3  are  now  eliminated  from  further  consideration, 
and  the  process  is  repeated  for  the  three  remaining  jobs 
and  machines  as  shown  in  Figure  10b.  Machine  1  now  has 
the  largest  difference,  and  it  is  assigned  to  job  1.  They 
are  also  eliminated  from  further  consideration  as  shown  in 
Figure  10c.  Machine  2  is  assigned  to  job  4  in  this  sequence. 
This  leaves  machine  3  to  be  assigned  to  job  2  by  default. 
These  job  assignments  are  optimal  in  that  the  cost  is  mini¬ 
mized  at  a  value  of  13. 

As  can  be  seen  in  the  example  problem,  this  method 
makes  assignments  on  the  basis  of  the  greatest  marginal 
cost  improvement,  i.e.,  it  takes  the  "greedy"  approach. 

In  the  example  problem  of  Figure  10a,  failure  to  assign 
machine  4  to  job  3  as  the  first  step  could  result  in  a  sub¬ 
sequent  cost  increase  of  3.  Rather  than  take  this  chance, 
the  "greedy"  method  makes  this  assignment  first  and  con¬ 
tinues  in  this  manner  until  all  assignments  are  made. 

Although  the  solution  was  optimal  for  this  example, 
this  method  does  not  guarantee  an  optimal  solution.  This 
occurs  for  two  reasons.  One  is  the  fact  that  ties  for  the 
largest  difference  are  broken  arbitrarily.  The  second, 
and  most  import  amt,  concerns  the  elimination  step.  As 
machines  and  jobs  are  assigned,  they  are  eliminated  from 
further  consideration.  This  precludes  their  use  in  any 
subsequent  tradeoffs  to  achieve  optimality. 
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The  "greedy"  method  uses  the  process  just  described 
to  assign  tankers  to  refueling  tracks  and  post -re fueling 
bases.  The  complete  methodology  is  illustrated  in 
Figure  11  and  discussed  below. 

The  first  step,  like  that  of  the  network  method, 
determines  how  many  tankers  should  be  assigned  to  each 
bomber.  The  next  step  computes  the  cost  of  assigning  every 
tanker  to  each  of  the  refueling  locations.  Like  the  net¬ 
work  method,  the  initial  refueling  locations  are  assumed 
to  be  those  of  an  average  tanker.  The  cost  in  this  case 
is  the  sum  of  the  fuel  consumed  to  reach  the  start  refuel¬ 
ing  point,  the  fuel  available  at  that  point,  and  the  fuel 
required  to  reach  the  post-refueling  base.  This  cost  is  com¬ 
puted  for  every  possible  tanker,  refueling  track,  and  PRB 
combination.  If  a  bomber  and  tanker  mating  is  infeasible 
due  to  timing,  the  costs  are  set  at  a  large  positive  value. 
After  the  costs  have  been  computed,  a  tanker  is  assigned  to 
a  particular  refueling  track  and  PRB  by  the  "greedy"  method. 
These  refueling  assignments  are  then  adjusted  to  optimize 
the  refueling  as  in  the  previous  methods.  The  final  step 
repeats  this  process  using  the  new  refueling  locations  as 
a  starting  point  and  checking  for  improvements. 

Computer 1 za t ion 

All  three  methods  were  programmed  in  FORTRAN  V  and 
run  on  a  Control  Data  Corporation  Cyber  750  computer. 
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Sample  inputs  and  outputs  from  these  programs  are  included 
in  Chapter  IV.  In  addition,  the  program  listings  and  sub¬ 
routines  are  included  in  Appendix  C.  Each  program  is  docu¬ 
mented  and  explained  by  means  of  comment  cards  which  are 
contained  in  the  program  listings. 

Verification  and  Validation 

The  validity  of  each  program  and  its  underlying 
methodology  were  evaluated  through  a  three-step  process 
adapted  from  the  work  of  Fishman  and  Kiviat  (Ref  5) .  The 
three  steps  of  this  process  are: 

1.  Verification  that  the  programs  work  as  designed 

2.  Validation  of  the  programs  against  real  world 
problems 

3.  Analysis  of  the  results 

The  verification  and  validation  steps  are  discussed  below. 
The  analysis  step  is  contained  in  Chapter  IV. 

Verification.  Proper  program  operation  was  verified 
by  insuring  that  the  main  functions  of  each  program  operated 
as  designed.  These  functions  included  flight  planning, 
determining  how  many  tankers  are  required  by  each  bomber, 
and  the  assignment  process.  Each  of  these  functions  are 
discussed  in  turn. 

The  flight  planning  function  was  checked  by  com¬ 
paring  computer  derived  figures  against  the  same  figures 
as  derived  from  the  performance  manuals.  This  was  done  for. 


all  computations  in  the  prototype  problem  and  for  selected 
computations  from  the  remaining  problems.  In  no  case  did 
the  computer  derived  figures  deviate  by  more  than  2  percent 
from  the  performance  manual  figures. 

The  function  of  determining  how  many  tankers  to 
assign  to  each  bomber  was  investigated  by  using  numerous 
example  problems  as  well  as  the  prototype  problem.  There 
was  no  occasion  where  this  function  did  not  assign  the 
proper  number  of  tankers  to  each  bomber. 

The  last  function  to  be  investigated  was  the  assign¬ 
ment  process.  Although  both  the  network  and  "greedy" 
methods  are  based  on  proven  algorithms  (GNET  and  "greedy"), 
they  were  still  checked  manually  against  example  problems 
and  the  prototype  problem  to  insure  that  they  worked  as 
designed.  This  proved  to  be  the  case. 

After  these  functions  were  verified  individually, 
the  complete  programs  were  verified  against  the  prototype 
problem  by  manually  checking  the  final  results.  The  final 
results  were  also  compared  against  each  other.  One  final 
indication  that  each  program  worked  as  designed  is  that  in 
those  cases  where. each  program  assigned  a  tanker  to  the 
same  refueling  track  and  PRB,  all  geographical  coordinates 
and  fuel  figures  were  identical.  If  this  had  not  been  the 
case,  it  would  have  indicated  a  fault  in  one  or  more  of 
the  programs. 
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Validation.  The  true  test  of  validity  for  any 
problem  solving  method  is  whether  or  not  it  can  solve  a 
real  world  problem.  This  test  was  not  applied  to  the 
methods  developed  for  this  study  because  of  the  actual 
problem's  high  degree  of  classification;  however,  these 
methods  were  used  to  solve  example  problems  that  were  care¬ 
fully  formulated  to  resemble  the  actual  problem.  They 
worked  as  designed  and  expected  against  these  problems  and 
demonstrated  face  validity  in  that  the  results  obtained 
were  entirely  reasonable.  Thus,  the  methods  developed  for 
this  study  were  validated  to  the  extent  that  the  sample 
problems  captured  the  real  world. 

Summary 

Three  methods  have  been  developed  to  solve  the 
bomber  and  tanker  mating  problem.  They  are  the  network 
method,  the  "greedy"  method,  and  the  pseudo-Logicon  method. 
These  methods  are  evaluated  against  several  problems  in  the 
next  chapter. 
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Chapter  IV.  Results  and  Analysis 

The  network,  "greedy,"  and  pseudo-Logicon  methods 
were  used  to  solve  five  bomber  and  tanker  mating  problems. 
These  problems  are  summarized  in  Table  VI  and  listed  in 
detail  in  Appendix  B.  The  increasing  sizes  of  these  prob¬ 
lems  reflect  the  "enrichment  and  elaboration"  process  that 
was  followed  in  the  development  stage  of  this  study.  All 
three  methods  were  developed  and  proven  against  the  proto¬ 
type  problem.  They  were  then  expanded  to  handle  the  larger 
problems  on  a  problem-by-problem  basis.  This  approach 
facilitated  the  programming,  debugging,  and  validation  of 
each  method.  The  final  goal  of  this  process  was  to  solve 
a  problem  the  size  of  the  normal,  day-to-day  alert  problem. 
This  goal  was  realized,  and  the  results  obtained  from  each 
method  are  reported  and  analyzed  in  this  chapter.  Problem 4 
served  as  the  primary  basis  of  comparison  because  it  repre¬ 
sented  the  alert  problem;  however,  the  results  from  the 
other  problems  were  studied  to  determine  if  one  method 
consistently  outperformed  the  others. 

Input  Data 

The  input  data  for  each  of  the  problems  consisted 
of  the  following  parameters: 
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1.  Geographical  coordinates  of  the  bomber,  tanker, 
and  PRBs 

2.  The  number  of  bombers  and  tankers  at  each  base 

3 .  The  tanker  level-off  gross  weights  for  each  base 

4 .  The  number  of  PRBs  and  their  capacities 

5.  Bomber  entry  point  fuel  requirements 

Some  of  this  data  is  included  in  Table  VI.  The  complete 
listing  for  each  problem  is  contained  in  Appendix  B . 


TABLE  VI 
PROBLEM  SUMMARY 


Problem 

Berber 

Bases 

Nunber  of 

Tanker 

Bases 

Nunber  of 
Tankers 

PRB 

Bases 

Prototype 

4 

4 

5 

7 

3 

1 

7 

10 

10 

17 

5 

2 

12 

26 

17 

39 

14 

- 

12 

52 

21 

78 

14 

4 

13 

90 

32 

135 

18 

Although  numbers  and  locations  may  vary  somewhat 
from  the  actual  figures  to  avoid  classification  difficulties, 
all  problems  other  than  the  prototype  problem  have  been 
structured  to  reflect  the  real  world.  This  was  accomplished 
through  the  following  techniques: 

1.  Bomber  and  tanker  bases  were  dispersed. through¬ 
out  the  United  States  in  general  geographic  areas  that  cor¬ 
respond  to  actual  bases. 


2.  Tanker  basing  reflects  active  duty.  National 
Guard,  and  Reserve  alert  force  commitments. 

3.  Tanker  gross  weights  vary  according  to  per¬ 
formance  limitations. 

4.  Bomber  to  tanker  ratios  correspond  to  actual 
figures  (Ref  5:72). 

5.  Post-refueling  bases  are  located  in  likely 
areas  such  as  Alaska,  Canada,  Greenland,  and  Iceland. 

Output  Data 

The  output  from  each  method  consists  of  two  parts. 
The  first  part  displays  each  bomber's  entry  point  fuel, 
the  deviation  from  required  entry  point  fuel,  and  the  total 
entry  point  fuel  for  all  bombers.  The  second  part  of  the 
output  lists  the  bomber  and  tanker  matings,  PRB  assignments, 
air  refueling  coordinates,  onload  and  offload  capabilities, 
and  time  on  the  refueling  track.  Sample  outputs  from  the 
network  method  solution  to  the  prototype  problem  are  dis¬ 
played  in  Figures  12  and  13.  Note  that  the  onload  and  off¬ 
load  capabilities  in  Figure  13  are  within  the  400-pound 
tolerance  established  in  Chapter  III. 

In  addition,  output  from  Subroutine  GNET  is  avail¬ 
able  if  desired.  This  output  data  includes  the  actual  arc 
flows  and  costs.  These  outputs  are  not  recommended  for 
larger  problems  because  they  quickly  become  voluminous. 
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Figure  13.  Sample  Output  from  Network  Method  and  Prototype  Problem,  Part 


Results 


The  results  of  running  the  network,  "greedy,"  and 
pseudo-Logicon  methods  against  each  of  the  five  problems 
are  summarized  in  Tables  VII  and  VIII.  Table  VII  displays 
the  results  obtained  from  the  first  iteration  for  each 
problem.  Table  VIII  displays  the  best  results  that  were 
obtained,  and  the  iteration  on  which  they  were  obtained. 

Best  in  this  case  is  defined  as  the  maximum  total  entry 
point  fuel. 

The  figures  in  these  tables  break  down  by  method 
the  number  of  bombers  that  arrive  at  the  entry  point  short 
of  the  required  fuel,  the  total  fuel  shortage,  and  the 
average  shortage  per  bomber.  Also  shown  are  the  number 
of  bombers  that  meet  or  exceed  the  required  entry  point 
fuel,  the  total  fuel  overage,  and  the  average  overage  per 
bomber.  The  final  column  displays  the  total  entry  point 
fuel  for  all  bombers. 

These  figures  were  selected  because  any  one  of 
them  can  be  used  as  evaluation  criteria  for  the  methods 
under  investigation;  however,  the  stated  objective  of  this 
study  was  to  develop  a  methodology  to  reduce  the  number  of 
bombers  requiring  degraded  tactics  and/or  to  reduce  the  dura¬ 
tion  of  these  tactics.  The  evaluation  criteria  that  cor¬ 
respond  to  this  goal  are  the  number  of  bombers  short  and 
the  average  shortage  per  bomber.  Thus,  the  method  that 
minimizes  both  of  these  criteria  will  obviously  be  the 
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FIRST  ITERATION  RESULTS 
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figures  in  1000' s  of  pounds 


TABLE  VZZ1 


figures  in  1000' a  of  pounds 


preferred  method.  If  both  criteria  are  not  minimized  by 
the  same  method,  then  the  method  that  also  maximizes  total 
entry  point  fuel  would  appear  to  have  the  advantage. 

A  review  of  Tables  VII  and  VIII  shows  that  against 
the  primary  problem  of  interest,  problem  4,  the  network 
method,  satisfies  all  three  criteria  for  the  single  and 
best  iteration  cases.  It  minimizes  the  number  of  shortages 
and  average  shortage  per  bomber  and  maximizes  the  total 
entry  point  fuel.  As  an  additional  check  for  consistency, 
it  satisfies  two  of  the  three  criteria  for  problems  1 
through  3  on  the  first  iteration  and  at  l°ast  two  of  the 
three  criteria  for  problems  2  and  3  on  the  best  iteration. 
The  only  other  method  to  satisfy  two  of  the  three  criteria 
is  the  "greedy"  method  on  the  prototype  problem  and 
problem  1.  These  were  the  only  inconsistent  results  noted 
and  are  most  likely  attributable  to  the  small  scale  of  the 
problems.  These  comparisons  are  summarized  in  Figure  14a. 

A  similar  comparison  of  the  "greedy"  and  pseudo- 
Logicon  methods  only  in  Figure  14b  shows  that  the  "greedy” 
method  satisfies  two  of  the  three  criteria  for  problem  4 
and  all  of  the  criteria  for  problem  3  in  both  cases.  For 
the  remaining  problems,  "greedy"  satisfies  a  minimum  of  two 
out  of  the  three  criteria. 

Based  on  these  comparisons,  the  network  method 
appears  to  be  the  best  method  for  solving  the  bomber  and 
tanker  mating  problem.  For  problem  4  it  reduces  the  number 
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First  Iteration  Best  Iteration 


Problem 

Bombers 

Short 

Average 

Shortage 

Maximum 
EP  Fuel 

Bombers 

Short 

Average 

Shortage 

Maximum 
EP  Fuel 

Prototype 

N/G/L 

G 

G 

N/G/L 

G 

G 

1 

N/G/L 

N 

G 

N/L 

G 

G 

2 

G 

N 

N 

N 

N 

N 

3 

N 

G 

N 

N 

G 

N 

4 

N 

N 

N 

N 

N 

N 

a.  Comparison  of  all  Three  Methods 


Problem 

First  Iteration 

Best  Iteration 

Bombers 

Short 

Average 

Shortage 

Maximum 
EP  Fuel 

Bombers 

Short 

Average 

Shortage 

Maximum 
EP  Fuel 

Prototype 

G/L 

G 

G 

G/L 

G 

G 

1 

G/L 

L 

G 

L 

G 

G 

2 

G 

L 

G 

G 

G 

G 

3 

G 

G 

G 

G 

G 

G 

4 

G 

L 

G 

G 

L 

G 

b.  Comparison  of  "Greedy"  and  Pseudo-Logicon  Methods 


NOTE :  Method  giving  the  best  results  for  each  criteria  where 

N  ■  Network  Method 
G  ■  "Greedy"  Method 
L  ”  Pseudo-Logicon  Method 


Figure  14.  Comparison  of  Methods 


of  bombers  shorted  by  17  percent  and  increased  total  entry 
point  fuel  by  16  percent  over  the  corresponding  figures 
for  the  pseudo-Log icon  method.  The  "greedy"  method  appears 
to  be  the  next  best  method  although  not  by  the  same  margin. 

It  reduces  the  number  of  bombers  shorted  by  7  percent  and 
increases  total  entry  point  fuel  by  8  percent  over  the 
corresponding  figures  for  the  pseudo-Log icon  method.  No 
method  showed  significant  advantages  in  reducing  the  average 
fuel  shortage  per  bomber. 

Analysis 

The  "greedy"  and  network  methods  appear  to  out¬ 
perform  the  pseudo-Logicon  method  because  the  latter  starts 
out  with  an  infeasible  solution,  i.e.,  it  ignores  the  PRB 
constraints.  It  then  has  to  go  back  and  send  the  excess 
tankers  to  unsaturated  bases  that  may  be  considerable  dis¬ 
tances  away.  The  end  result  is  that  these  tankers'  off¬ 
load  capabilities  are  adversely  affected  which  in  turn 
adversely  affects  the  bombers'  entry  point  fuel. 

The  "greedy"  method,  unlike  the  pseudo-Logicon 
method,  deals  only  with  feasible  solutions;  however,  as  pre¬ 
viously  noted,  once  it  makes  a  tanker,  bomber,  and  PRB 
assignment  it  is  unable  to  go  back  and  perform  the  necessary, 
tradeoffs  to  improve  the  solution. 


The  network  method  works  best  because  it  takes  all 
of  the  constraints  (including  PRB  capacities)  into  account. 


deals  only  with  feasible  solutions,  and  can  perform  the 
necessary  tradeoffs  through  the  primal-simplex  method  to 
improve  the  solution.  This  is  reflected  in  the  results 
that  were  obtained. 

Effects  of  Iteration 

The  prototype  problem  was  so  small  that  iterating 
had  no  effect  on  any  of  the  methods  (see  Tables  VII  and 
VIII) .  Each  iteration  resulted  in  the  same  bomber,  tanker, 
and  PRB  assignments.  Iteration  of  the  larger  problems  did 
result  in  reassignment  of  some  tankers  to  different  bombers 
or  PRBs.  This  was  an  expected  result  because  of  the  changed 
air  refueling  locations.  The  result  that  was  not  expected 
was  the  small  improvements,  if  any,  in  the  evaluation  cri¬ 
teria.  In  some  cases,  all  three  criteria  were  improved. 

The  "greedy"  solution  for  problem  4  is  one  such  example. 

In  other  cases  no  improvements  were  noted  such  as  the  net¬ 
work  solution  for  problem  4 .  There  were  also  cases  where 
the  results  were  mixed  as  in  the  "greedy"  solution  for 
problem  2.  Finally,  there  were  some  cases  where  all  three 
criteria  actually  decreased;  however,  this  is  not  shown  in 
Table  VIII  since  it  reflects  the  best  iteration. 

Two  factors  appear  to  be  responsible  for  these 
inconsistent  results.  One  is  the  large  variation  in  tanker 
gross  weights.  These  variations  range  from  250,000  pounds 
to  279,500  pounds  at  level -off.  The  other  factor  involved 
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is  the  variable  distances  to  the  PRBs.  These  two  factors 
combine  to  produce  a  large  number  of  tankers  that  deviate 
significantly  from  the  average  tankers.  These  variations 
in  turn,  can  require  large  shifts  in  air  refueling  loca¬ 
tions  in  order  to  optimize  individual  air  refuelings. 

These  large  shifts  can  significantly  alter  the  arc  costs 
from  iteration  to  iteration,  and  there  is  no  guarantee  that 
the  minimum  cost  flow  on  a  subsequent  iteration  will  be  less 
than  the  minimum  cost  on  the  current  iteration.  Some  refuel¬ 
ing  tracks  are  moved  away  from  the  bomber  and  tanker  bases, 
and  some  are  moved  in  the  opposite  direction.  If  the  refuel¬ 
ing  track  is  moved  away  from  the  tanker's  departure  base, 
the  tanker  cost  to  get  to  the  ARCP  is  increased.  In  addi¬ 
tion,  the  tanker's  cost  from  the  EAR  point  to  the  PRB  may 
also  increase  due  to  this  adjustment.  This  increases  the 
total  cost  for  that  tanker.  With  the  large  numbers  of 
bombers  and  tankers  involved,  the  next  iteration  may  have 
a  higher  total  cost.  The  net  result  of  this  argument  is 
that  contrary  to  our  initial  hypothesis,  minimizing  the 
cost  of  refueling  does  not  necessarily  produce  the  best  entry 
point  solution,  but  does  produce  a  good  one. 

Effect  of  Tanker  Constraints 

The  tanker  inputs  for  eac*.  method  include  the  base 
location,  number  of  aircraft  at  a  base,  and  gross  weight. 


Each  of  these  factors  can  effect  the  mating  procedure  to  a 
certain  extent. 

Location  and  Number.  These  factors  do  not  affect 
the  problem  significantly  because  all  aircraft  are  assumed 
to  take-off  at  the  same  time.  This  also  means  that  they 
are  airborne  for  the  same  amount  of  time.  Thus,  the  only 
advantage  gained  from  being  closer  to  the  ARCP  is  that  an 
aircraft  may  be  able  to  hold  at  maximum  endurance  airspeed 
while  waiting  for  a  bomber.  Holding  saves  fuel  but  the 
difference  is  not  significant.  For  example,  assume  that 
two  identical  tankers  are  assigned  to  ARCPs  1000  and  2000 
miles  from  takeoff  with  start  air  refueling  times  corres¬ 
ponding  to  the  time  it  takes  to  fly  the  2000  miles.  The 
tanker  that  has  to  fly  directly  to  the  ARCP  consumes 
111,500  pounds  of  fuel.  The  tanker  that  flies  1000  miles 
and  holds  until  the  first  tanker  reaches  its  ARCP  will 
consume  110,200  pounds  of  fuel.  Thus,  the  advantage  gained 
from  being  1000  miles  closer  to  the  ARCP  is  only  1300  pounds. 

Gross  Weight.  As  previously  discussed,  large  vari¬ 
ations  in  the  tanker  gross  weights  result  in  large  shifts 
in  the  air  refueling  locations  and  costs.  This  factor  was 
investigated  by  making  all  tanker  gross  weights  equal. 

When  this  was  done,  the  network  method  still  obtained  the 
best  results  followed  by  the  "greedy"  method. 
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Effect  of  PRB  Constraints 

The  PRB  locations  and  capacities  were  the  critical 
factors  in  this  investigation.  Routes  to  the  entry  points 
and  locations  of  the  PRBs  caused  certain  PRBs  to  be  favored 
over  others.  For  example,  out  of  18  possible  PRBs  in 
problem  4 ,  only  6  were  used  when  all  capacity  constraints 
were  removed.  One  of  these  had  72  aircraft  assigned  when 
its  capacity  was  only  11.  This  demonstrates  why  the  net¬ 
work  and  "greedy"  methods  obtain  better  solutions  than  the 
pseudo-Log icon  method.  It  has  to  move  61  aircraft  to  new 
PRBs  at  an  obviously  large  penalty. 

Summary 

The  network,  "greedy,"  and  pseudo-Logicon  methods 
were  evaluated  against  five  mating  problems  in  this  chapter. 
Input  and  output  data  were  described  and  the  results  were 
reported  and  analyzed.  The  resulting  conclusions  and  recom¬ 
mendations  are  presented  in  the  next  chapter. 


Chapter  V.  Conclusions  and  Recommendations 

Summary 

The  purpose  of  this  research  effort  was  to  investi¬ 
gate  the  current  methodology  used  to  mate  bombers  and 
tankers  in  the  Single  Integrated  Operations  Plan  with  an 
objective  of  improving  the  process  if  possible.  Two  methods 
were  formulated  to  achieve  this  objective.  One  used  net¬ 
work  theory  in  an  attempt  to  obtain  an  optimal  solution. 

The  second  used  a  "greedy"  method  to  provide  a  feasible  but 
not  necessarily  optimal  solution  as  an  alternative  approach. 
Both  of  these  methods  were  then  compared  to  the  method  cur¬ 
rently  in  use.  This  comparison  was  based  on  five  problems 
of  progressively  increasing  difficulty,  concluding  with  a 
problem  that  was  structured  to  reflect  an  actual  SIOP 
mating  problem. 

Conclusions 

The  size  and  complexity  of  the  bomber  and  tanker 
mating  problem  precludes  a  truly  optimal  solution.  The 
interdependence  of  bomber  and  tanker  assignments,  air  refuel 
ing  locations,  and  PRB  assignments  result  in  a  problem  com¬ 
monly  referred  to  as  a  three-dimensional  assignment  problem. 
There  is  no  known  polynomial  bounded  algorithm  for  solving 
such  a  problem. 


The  network  method  was  reformulated  to  fix  the  air 


refueling  locations,  obtain  an  initial  feasible  solution, 
and  then  iterate  this  solution  to  improve  it.  This  revised 
method  proved  to  be  the  best  of  the  three  methods  under 
investigation  followed  by  the  "greedy"  method. 

Iterating  the  three  methods  did  not  result  in  any 
significant  improvements  to  the  initial  solutions. 

As  long  as  PRB  locations  and  capacities  are  not  a 
factor,  the  current  method  is  essentially  identical  to  the 
network  method;  however,  when  these  constraints  are  a  fac¬ 
tor,  the  current  method  is  penalized  because  it  does  not 
take  them  into  account  until  after  the  initial  assignments 
have  been  made. 

Recommendations 

The  Strategic  Air  Command  should  investigate  the 
possibility  of  incorporating  PRB  capacities  in  the  network 
solving  algorithm  employed  in  the  Mating  and  Ranging  Pro¬ 
gram. 

Recommended  Areas  for  Follow-on  Study 


The  network  method  developed  in  this  study  should 
be  expanded  to  include  some  sort  of  bomber  priority  in  being 
assigned  a  tanker.  In  this  way,  a  bomber  short  on  EP  fuel 
would  be  given  a  higher  priority  for  being  assigned  a 
strong  tanker.  Those  bombers  over  their  EP  fuel  require¬ 
ments  would  be  assigned  a  lower  priority. 


A  program  should  be  developed  to  relate  a  bomber's 
fuel  state  to  his  probability  of  survival  through  enemy 
territory.  For  example,  a  bomber  that  meets  or  exceeds 
its  EP  fuel  requirements  would  be  given  a  survival  proba¬ 
bility  of  1.0,  while  those  not  meeting  their  requirements 
would  have  some  lower  probability  that  would  depend  on  the 
bomber ' s  route  of  flight  and  the  additional  enemy  defenses 
encountered.  As  a  result  of  degraded  tactics,  such  a  pro¬ 
gram  could  then  assign  tankers  based  on  bomber  fleet  sur¬ 
vivability  rather  than  explicit  fuel  requirements. 

A  study  should  be  undertaken  to  determine  the 
effect  of  using  KC-lOs  and  re-engined  KC-135s  in  the  SIOP. 
These  aircraft  offer  greater  fuel  offload  capabilities 
than  the  KC-135A,  and  could  increase  bomber  entry  point 
fuel  substantially. 

Comments 

The  addition  of  cruise  missile  commitments  to  the 
B-52  fleet  in  the  early  1980s  has  a  three-fold  effect  on 
the  tanker  assignment  problems.  First,  the  addition  of 
cruise  missiles  on  the  aircraft  decreases  the  bomber's  fuel 
carrying  capability.  Secondly,  the  increased  drag  from 
these  missiles  increases  fuel  consumption.  Finally,  as  the 
cruise  missile  carrying  aircraft  assume  a  stand-off  role, 
they  will  most  likely  recover  into  bases  that  are  currently 
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used  as  tanker  PRBs.  This  means  that  both  tankers  and 
bombers  will  be  competing  for  the  PRB  space.  The  end 
result  of  these  effects  is  that  the  tanker  to  bomber  mating 
problem  will  become  even  more  critical. 
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The  general  model  used  in  deriving  the  aircraft 
fuel  flow  equations  is  the  following: 

FF  »  bQ  *  GW  +  bL 

where  bg  and  b^  are  constants,  FF  is  the  estimated  fuel 
flow  in  pounds  per  minute,  and  GW  is  the  aircraft  gross 
weight  in  thousands  of  pounds.  It  was  possible  to  model 
aircraft  fuel  flow  in  this  manner  since  it  was  assumed  that 
the  aircraft  fly  a  maximum  range  cruise  profile,  and  there¬ 
fore  their  fuel  flows  depend  only  on  gross  weight  changes. 
Endurance  fuel  flows  can  be  modeled  the  same  way  since  a 
constant  altitude  (30,000  feet)  is  assumed,  and  air  refuel¬ 
ing  fuel  flows  are  identical  to  cruise  fuel  flows  at  con¬ 
stant  altitude  and  airspeed  with  the  addition  of  a  fuel 
flow  degradation  factor.  The  end  result  of  all  of  these 
factors  is  that  all  fuel  flows  are  dependent  only  on  changes 
in  aircraft  gross  weight,  and  simple  linear  regressions 
can  be  performed  for  each  different  phase  of  flight  for 
each  aircraft. 

The  maximum  range  cruise  (MRC)  fuel  flows  for  both 
aircraft  were  divided  into  two  gross  weight  categories. 

The  B-52  used  those  weights  above  340,000  pounds  as  one 
category,  and  those  weights  equal  to  or  below  this  weight. 
The  KC-135  used  180,000  pounds  as  the  dividing  point.  The 
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two  weight  categories  were  used  to  provide  a  better  linear 
estimation  of  the  fuel  flows.  The  two  specific  weights 
selected  as  dividing  points  were  chosen  because  they 
represent  weights  typical  of  mean  values  encountered 
throughout  the  mission  profile. 

Table  A-l  summarizes  the  data  used  in  the  KC-135 
regressions,  and  Table  A- 3  contains  the  B-52  data.  All 
of  the  data  points  were  extracted  from  the  appropriate 
aircraft  performance  manual,  and  used  a  standard  tempera¬ 
ture  deviation  of  0 . 

The  Statistical  Package  for  the  Social  Sciences 
(SPSS)  linear  regression  routine  was  used  to  derive  all 
seven  fuel  flow  equations,  and  the  resulting  equations 
are  shown  in  Table  A-5. 

SPSS  summary  tables  for  all  equations  are  given 
in  Tables  A- 2  and  A-4 .  Statistically,  all  of  the  regres¬ 
sion  models  are  highly  significant,  with  the  lowest  coeffi- 

2 

cient  of  determination  (R  )  value  being  .988.  This  indi¬ 
cates  that  almost  99  percent  cf  the  variability  of  fuel 
flow  is  explained  by  the  regression  model.  The  high  over¬ 
all  F  values  obtained  in  all  cases  cor  firms  that  gross 
weight  contributes  significantly  to  the  regression  models. 
All  residuals  (the  difference  between  actual  and  predicted 
values)  were  within  two  standard  deviations  of  the  mean 
response,  again  indicating  the  validity  of  the  model. 
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TABLE  A-l 


KC-135  DATA  USED  FOR  FUEL  FLOW  REGRESSIONS 


ENDURANCE  AT  30,000  FEET 

GW  (1000s  of  lbs) 

FF  (lbs/min) 

260 

217.8 

240 

200.0 

220 

181.3 

200 

163.3 

180 

145.9 

MRC  FOR  GROSS  WEIGHTS  GREATER  THAN  180,000  POUNDS 
GW  (1000s  of  lbs)  FF  (lbs/min) 


280 

243.5 

270 

233.7 

260 

225.3 

250 

217.5 

240 

207.9 

230 

198.6 

22  0 

191.0 

210 

182.8 

200 

174.3 

190 

165.9 

180 

157.2 

TABLE  A-l — Continued 


MRC  FOR  GROSS  WEIGHTS  LESS  THAN  OR  EQUAL  TO 
180.000  POUNDS 


GW  (10008  of  lbs) 


FF  (lbs/min) 


180 

170 

160 

150 

140 

130 

120 


157.2 

149.2 

140.6 

132.7 
124.6 

115.8 
109.1 


AIR  REFUELING 


GW  (10003  of  lbs) 


FF  (lbs/min) 


210 

180.0 

200 

173.4 

190 

165.9 

180 

160.9 

170 

154.8 

160 

150.5 

150 

145.5 

140 

141.6 

130 

138.2 

TABLE  A- 2 


SPSS  SUMMARY  TABLES  FOR  THE  KC-135A 


ENDURANCE  AT  30,000  FEET 


Overall  F 
Significance 
Multiple  R 
R  Square 


22890.1 

.000 

.999 

.999 


I  MRC  FOR  GROSS  WEIGHTS  GREATER  THAN  180,000  POUNDS 


Overall  F 
Significance 
Multiple  R 
R  Square 


22662.5 

.000 

.999 

.999 


MRC  FOR  GROSS  WEIGHTS  LESS  THAN  OR  EQUAL  TO 
180,000  POUNDS 


Overall  F  6911.1 

Significance  .000 

Multiple  R  .999 

R  Square  .999 


AIR  REFUELING 


Overall  F  553.6 
Significance  .000 
Multiple  R  .993 
R  Square  .988 


TABLE  A-3 


B-52H  DATA  USED  FOR  FUEL  FLOW  REGRESSIONS 


MRC  FOR  GROSS  WEIGHTS  GREATER  THAN 
340,000  POUNDS 

GW  (1000s  of  lbs)  FF  (lbs /min) 


4  76 

350.3 

4  70 

348.5 

460 

341.6 

450 

333.8 

440 

326.0 

430 

317.3 

420 

312.1 

410 

305.2 

400 

296.5 

390 

289.6 

380 

283.5 

370 

277.4 

360 

268.6 

350 

261.8 

340 

253.2 

MRC  FOR  GROSS  WEIGHTS  LESS  THAN  OR  EQUAL  TO 
340,000  POUNDS 


GW  (1000s  of  lbs)  FF  (lbs/min) 


340 

253.2 

330 

248.8 

320 

242.8 

310 

235.8 

300 

229.8 

290 

220.0 
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TABLE  A-4 

SPSS  SUMMARY  TABLES  FOR  THE  B-52H 


MRC  FOR  GROSS  WEIGHTS  GREATER  THAN 
340.000  POUNDS 


Overall  F 
Significance 
Multiple  R 
R  Square 


13814.1 

.000 

.999 

.999 


MRC  FOR  GROSS  WEIGHTS  LESS  THAN  OR  EQUAL  TO 


340j_000  POUNDS 


Overall  F 
Significance 
Multiple  R 
R  Square 


667.3 

.000 

.997 

.994 


AIR  REFUELING 


Overall  F  2263.5 

Significance  .000 

Multiple  R  .998 

R  Square  .996 


TABLE  A- 5 

SUMMARY  OF  REGRESSION  EQUATIONS 


KC-135A: 

Endurance 

MRCC  GW  GT  180,000  lbs 
MRC0  GW  LE  180,000  lbs 
Air  Refueling 


FF».9025*GW-16.89 
FF*.8564*GW+  2.83 
FF«.8132*GW+10.80 
FF* . 5237*GW+67 . 71 


B-52H: 

MRC@  GW  GT  340,000  lbs 
MRC@  GW  LE  340,000  lbs 
Air  Refueling 


FF=*.  7178*GW+10 .27 
FF*.6286*GW+40.73 
FF=.7837*GW+90.53 


FF  ■  Fuel  Flow  in  lbs/min 

GW  *  Aircraft  Gross  Weight  in  1000s  of  lbs 


This  appendix  contains  the  data  sets  for  all  of  the 
scenarios  used  in  this  report.  Data  sets  are  listed  as 
they  were  input  to  each  of  the  models.  Bomber  input  data 
was  stored  on  tape  1  in  the  following  sequence:  departure 
latitude,  departure  longitude,  entry  point  latitude, 
entry  point  longitude,  and  entry  point  fuel  desired. 

Tanker  input  data  was  stored  on  tape  2  in  the  following 
sequence:  departure  latitude,  departure  longitude,  gross 
weight  at  level-off,  and  number  of  tankers  at  this  loca¬ 
tion.  Recovery  base  input  data  was  stored  on  tape  3  in 
this  sequence:  latitude,  longitude,  and  capacity.  All 
latitudes  and  longitudes  were  input  as  degrees  and  frac¬ 
tions  of  degrees  (e.g.,  35  30*  was  input  as  35.5).  Eastern 
longitudes  and  southern  latitudes  are  input  as  negative 
numbers  (e.g.,  35  E  is  input  as  -35.). 


PROTOTYPE  PROBLEM 


BOMBER  ORTA: 


Departure 

Departure 

Entry  Point 

Entry  Point 

Entry  Point 

Latitude 

Longitude 

Latitude 

Longitude 

Fuel  Desired 

35.00 

110.00 

85.00 

-150.00 

205 

18.00 

100.00 

87.00 

-160.00 

250 

30.00 

91.00 

88.00 

140.00 

220 

20.00 

80.00 

88.00 

50.00 

250 

TANKER  DATA: 


Departure  Departure 

Latitude  Longitude 


Level-Off 

Gross  Weight  Number 


35.00 

110.00 

37.00 

105.00 

18.00 

100.00 

30.00 

91.00 

35.00 

85.00 

279.5  1 

279.5  2 

279.5  1 

279.5  1 

264.5  2 


RECOVERY  BASE  DATA: 


Latitude 

65.00 

62.00 

65.00 


Longitude 

115.00 

100.00 

85.00 


Capacity 

2 

2 

10 
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PROBLEM  ONE 


BOMBER  DATA: 

Departure  Departure 

T-*titude  Longitude 


Entry  Point  Entry  Point  Entry  Point 

Latitude  Longitude  Puel  Desired 


38.00 

45.00 

43.00 

43.00 

33.00 

46.00 

46.00 

32.00 

32.00 

44.00 


121.00 

115.00 

105.00 

105.00 

100.00 

95.00 

95.00 

94.00 

94.00 

84.00 


78.31 

82.38 

79.57 

80.86 

76.88 

88.78 

86.98 

77.03 

83.16 

82.49 


-119.81  260 

175.50  230 

-98 . 72  213 

-90.82  220 

-88.31  255 

-49.62  233 

-56.22  227 

-64.60  253 

15.96  250 

18.29  233 


TANKER  DATA: 

Departure  Departure  Level-Off 

Latitude  Longitude  Gross  Weight 

38.00  121.00  279.5 

45.00  115.00  277.5 

43.00  105.00  259.5 

33.00  100.00  279.5 

46.00  95.00  279.5 

32.00  94.00  279.5 

44.00  84.00  279.5 

42.00  87.00  268.5 

43.00  87.00  279.5 

41.00  112.00  259.5 


Number 

2 

2 

2 

1 

2 

2 

2 

1 

1 

1 


RECOVERY  BARE  DATA: 

Latitude 

Longitude 

Capacity 

69.00 

50.00 

5 

65.00 

145.00 

5 

75.00 

55.00 

5 

55.00 

115.00 

3 

52.00 

107.00 

3 
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PROBLEM  TWO 

P 

\\ 

BOMBER  DATA: 

r 

Departure 

Departure 

Entry  Point 

Entry  Point 

!  * 

Entry  Point 

Latitude 

Longitude 

Latitude 

Longitude 

Fuel  Desired 

38-00 

121.00 

78.30 

-119.80 

260 

38.00 

121.00 

80.70 

-133.60 

220 

45.00 

115.00 

83.60 

-158.20 

240  !i 

45.00 

115.00 

82.60 

-145.40 

230  {I 

43.00 

105.00 

79.60 

-98.70 

255 

43.00 

105.00 

86.80 

-142.10 

222 

43.00 

105.00 

86.90 

148.80 

240 

47.00 

100.00 

88.10 

144.60 

235 

48.00 

100.00 

89.20 

-144.40 

230  :  ‘ 

46.00 

95.00 

88.80 

-49.60 

230 

46.00 

95.00 

85.40 

-65.90 

240  j 

46.00 

88.00 

83.60 

-70.90 

218  ! 

46.00 

88.00 

85.50 

-51.30 

225 

44.00 

84.00 

85.70 

-36.80 

220 

44.00 

84.00 

83.40 

-51.10 

230 

33.00 

100,00 

77.00 

-84.20 

255 

33.00 

100.00 

85.80 

-69.30 

245 

33.00 

97.00 

86.60 

-52.90 

245 

33.00 

97.00 

77.50 

-69.90 

255  i: 

32.00 

94.00 

77.70 

-77.10 

250  i 

32.00 

94.00 

86.20 

-42.90 

245  ! 

32.00 

94.00 

86.40 

24.20 

250  j: 

35.00 

90.00 

83.30 

83.60 

230 

35.00 

90.00 

85.50 

69.60 

225 

32.00 

85.00 

83.40 

24.10 

245 

32.00 

85.00 

81.50 

6.70 

220 
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problem  TWO — Continued 


TANKER  DATA: 

Departure 

Depature 

Latitude 

Longitude 

38.00 

121.00 

45.00 

115.00 

43.00 

105.00 

47.00 

100.00 

46.00 

95.00 

46.00 

88.00 

44.00 

84.00 

33.00 

100.00 

33.00 

97.00 

32.00 

94.00 

35.00 

90.00 

32.00 

85.00 

38.00 

122.00 

41.00 

112.00 

43.00 

87.00 

42.00 

87.00 

36.00 

84.00 

RECOVERY  base 

DATA: 

Latitude 

Longitude 

62.00 

150.00 

65.00 

145.00 

55.00 

115.00 

52.00 

107.00 

75.00 

55.00 

69.00 

50.00 

65.00 

155.00 

65.00 

157.00 

49.00 

54.00 

52.00 

60.00 

64.00 

68.00 

56.00 

111.00 

54.00 

110.00 

65.00 

20.00 

Level-Off 

Gross  Weight  Number 


279.5  3 

277.5  3 

259.5  2 

279.5  3 

279.5  3 

276.5  3 

279.5  2 

279.5  2 

276.3  3 

279.5  3 

279.5  3 

279.5  2 

278.5  3 

259.5  1 

279.5  1 

268.5  1 

250.0  1 


Capacity 

5 

3 

3 

3 

3 

3 

1 

2 

3 

3 

3 

3 

3 

5 
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PROBLEM  THREE 


BOMBER  DATA: 

Departure 

Latitude 

38.00 

38.00 

38.00 

38.00 

45.00 

45.00 

45.00 

45.00 

43.00 

43.00 

43.00 

43.00 

43.00 

43.00 

47.00 

47.00 

47.00 

47.00 

46.00 

46.00 

46.00 

46.00 

46.00 

46.00 

46.00 

46.00 

44.00 

44.00 

44.00 

44.00 

33.00 

33.00 

33.00 

33.00 

33.00 

33.00 

33.00 

33.00 

32.00 

32.00 

32.00 

32.00 


Departure 

Longitude 

121.00 

121.00 

121.00 

121.00 

115.00 

115.00 

115.00 

115.00 

105.00 

105.00 

105.00 

105.00 

105.00 

105.00 

100.00 

100.00 

100.00 

100.00 

95.00 

95.00 

95.00 

95.00 

88.00 

88.00 

88.00 

88.00 

84.00 

84.00 

84.00 

84.00 

100.00 

100.00 

100.00 

100.00 

97.00 

97.00 

97.00 

97.00 

94.00 

94.00 

94.00 

94.00 


Entry  Point 
Latitude 

78.30 
79.00 

80.70 

80.60 

83.60 

82.60 
83.00 

83.70 

83.30 

84.90 
85.00 

79.60 

86.80 

86.90 

85.90 

88.10 

87.70 

89.20 

88.80 

85.40 

83.60 

85.50 

85.40 

83.60 

85.50 
87.00 

85.70 

84.30 

83.40 

82.90 

77.70 
77.00 

85.30 

86.70 
77.00 

86.60 

87.50 

77.50 

77.70 

86.20 

78.40 

86.40 


Entry  Point 
Longitude 

-119.80 

-155.10 

-133.60 

-122.10 

-158.20 

-145.40 

-129.60 

-134.40 

-141.10 

-171.10 

141.10 

-98.70 

-142.10 

148.80 

-106.60 

144.60 

-108.10 

-144.40 

-49.60 

-65.90 

-63.90 

-44.30 

-72.90 

-70.90 

-51.30 

-63.20 

-36.80 

-48.50 

-51.10 

-46.30 

-88.80 

-84.20 

-69.30 

-56.50 

-78.80 

-52.90 

-24.60 

-69.90 

-77.10 

-42.90 

-69.50 

24.20 


Entry  Point 
Fuel  Desired 

260 

220 

220 

250 

240 

230 

230 

224 
230 

225 
220 
255 
222 
240 
235 
235 
230 
230 
230 
240 
240 
220 
220 
218 
225 
220 
220 
228 
230 
225 
260 

255 
245 

249 
258 
245 

256 
255 

250 
245 

251 
250 


96 


A 


BOMBER  DATA— -Continued 


Departure 

Departure 

Entry  Point 

Entry  Point 

Entry  Point 

Latitude 

Longitude 

Latitude 

Longitude 

Fuel  Desired 

32.00 

94.00 

86.20 

.20 

250 

32.00 

94.00 

85.30 

-19.70 

240 

35.00 

90.00 

83.30 

83.60 

230 

35.00 

90.00 

85.50 

69.60 

225 

35.00 

90.00 

84.60 

69.10 

225 

35.00 

90.00 

83.70 

68.80 

220 

32.00 

85.00 

83.40 

24.10 

245 

32.00 

85.00 

80.10 

-47.90 

250 

32.00 

85.00 

82.40 

-5.10 

250 

32.00 

85.00 

81.50 

6.70 

220 

TANKER  DATA: 


iparture 

Departure 

Level-Off 

ititude 

Longitude 

Gross  Height 

Number 

38.00 

121.00 

279.5 

4 

45.00 

115.00 

277.5 

5 

43.00 

105.00 

259.5 

5 

47.00 

100.00 

279.5 

5 

46.00 

95.00 

279.5 

5 

46.00 

88.00 

276.5 

5 

44.00 

84.00 

279.5 

5 

33.00 

100.00 

279.5 

4 

33.00 

97.00 

276.3 

4 

32.00 

94.00 

279.5 

4 

35.00 

90.00 

279.5 

4 

32.00 

85.00 

279.5 

4 

38.00 

122.00 

278.5 

4 

41.00 

112.00 

259.5 

1 

43.00 

87.00 

279.5 

1 

42.00 

87.00 

268.5 

1 

36.00 

84.00 

250.0 

1 

41.00 

86.00 

279.5 

5 

38.00 

120.00 

268.0 

4 

35.00 

99.00 

279.5 

4 

37.00 

120.00 

279.5 

3 
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BO*BER  DATA: 


Departure 

Departure 

Entry  Point 

Entry  Point 

Entry  Point 

Latitude 

Longitude 

Latitude 

Longitude 

Fuel  Desired 

38.00 

121.00 

78.10 

-125.50 

250 

38.00 

121.00 

79.30 

-119.80 

260 

38.00 

121.00 

79.00 

-155.10 

220 

38.00 

121.00 

80.70 

-133.60 

220 

38.00 

121.00 

80.60 

-122.10 

250 

45.00 

115.00 

85.00 

-148.00 

225 

45.00 

115.00 

83.60 

-158 . 20 

240 

45.00 

115.00 

82.60 

-145.40 

230 

45.00 

115.00 

83.00 

-129.60 

230 

45.00 

115.00 

83.70 

-134.40 

224 

43.00 

105.00 

83.20 

-156.70 

228 

43.00 

105.00 

83.70 

-163.30 

224 

43.00 

105.00 

83.30 

-141.10 

230 

43.00 

105.00 

84.90 

-171.90 

225 

43.00 

105.00 

85.00 

141.10 

220 

43.00 

105.00 

79.60 

-98.70 

255 

43.00 

105.00 

86.80 

-142.10 

222 

43.00 

105.00 

86.90 

148.80 

240 

47.00 

100.00 

85.40 

-70.70 

235 

47.00 

100.00 

85.90 

-106.60 

235 

47.00 

100.00 

88.10 

144.60 

235 

47.00 

100.00 

87.70 

-108.10 

230 

47.00 

100.00 

89.20 

-144.40 

230 

46.00 

95.00 

84.40 

-61.10 

220 

46.00 

95.00 

88.80 

-49.60 

230 

46.00 

95.00 

85.40 

-65.90 

240 

46.00 

95.00 

83.60 

-63.90 

240 

46.00 

95.00 

85.50 

-44.30 

220 

46.00 

88.00 

84.80 

-83.60 

225 

46.00 

88.00 

85.40 

-72.90 

220 

46.00 

88.00 

83.60 

-70.90 

218 

46.00 

88.00 

85.50 

-51.30 

225 

46.00 

88.00 

87.00 

-63.20 

220 

44.00 

84.00 

87.40 

39.40 

220 

44.00 

84.00 

85.70 

-36.80 

220 

44.00 

84.00 

84.30 

-48.50 

228 

44.00 

84.00 

83.40 

-51.10 

230 

44.00 

84.00 

82.90 

-46.30 

225 

37.00 

120.00 

82.90 

-128.80 

251 

37.00 

120.00 

84.00 

-132.60 

252 

37.00 

120.00 

83.90 

-115.00 

250 

37.00 

120.00 

86.00 

-146.70 

220 

34.00 

117.00 

81.60 

-151.20 

251 

34.00 

117.00 

81.60 

173.60 

220 

34.00 

117.00 

84.20 

-164.80 

255 
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BOMBER  DATA- -Continued 


Departure  Departure 

Latitude  Longitude 


Entry  Point 
Latitude 


Entry  Point 
Longitude 


Entry  Point 
Fuel  Desired 


34.00 

117.00 

84.40 

176.70 

222 

34.00 

117.00 

85.50 

-178.50 

257 

33.00 

100.00 

75.90 

-91.70 

250 

33.00 

100.00 

88.20 

-123.90 

256 

33.00 

100.00 

77.70 

-88.80 

260 

33.00 

100.00 

77.00 

-84.20 

255 

33.00 

100.00 

85.30 

-69.30 

245 

33.00 

100.00 

86.70 

-56.50 

249 

33.00 

97.00 

88.70 

22.30 

252 

33.00 

97.00 

86.60 

-.40 

250 

33.00 

97.00 

83.90 

-42.20 

250 

33.00 

97.00 

84.30 

-33.50 

240 

33.00 

97.00 

77.00 

-78.80 

258 

33.00 

97.00 

86.60 

-52.90 

245 

33.00 

97.00 

87.50 

-24.60 

256 

33.00 

97.00 

77.50 

-69.90 

255 

32.00 

94.00 

78.00 

-90.50 

252 

32.00 

94.00 

84.60 

-76.50 

251 

32.00 

94.00 

83.70 

-29.70 

254 

32.00 

94.00 

77.70 

-77.10 

250 

32.00 

94.00 

86.20 

-42.90 

245 

32.00 

94.00 

78.40 

-69.50 

251 

32.00 

94.00 

86.40 

24.20 

250 

32.00 

94.00 

86.20 

.20 

250 

32.00 

94.00 

85.30 

-19.70 

240 

35.00 

90.00 

86.60 

103 . 20 

219 

35.00 

90.00 

83.30 

83.60 

230 

35.00 

90.00 

85.50 

69.60 

225 

35.00 

90.00 

84.60 

69.10 

22S 

35.00 

90.00 

83.70 

68.80 

220 

32.00 

85.00 

84.90 

2.20 

250 

32.00 

85.00 

83.40 

24.10 

245 

32.00 

85.00 

80.10 

-47.90 

250 

32.00 

85.00 

82.40 

-5.10 

250 

32.00 

85.00 

81.50 

6.70 

220 

43.00 

75.00 

89.20 

9.90 

225 

43.00 

75.00 

85.90 

-82.80 

220 

43.00 

75.00 

88.10 

-29.40 

230 

43.00 

75.00 

86.30 

-67.20 

218 

43.00 

75.00 

86.80 

-37.90 

220 

46.00 

68.00 

87.20 

-96.70 

223 

46.00 

68.00 

87.70 

-73.70 

227 

46.00 

68.00 

85.20 

-84.40 

230 

46.00 

68.00 

84.10 

-81.30 

235 

46.00 

68.00 

85.60 

-57.10 

225 

100 


TANKER  DATA 


Departure 

Departure 

Level-Off 

Latitude 

Longitude 

Gross  Weight 

Number 

38.00 

121.00 

279.5 

5 

37.00 

120.00 

279.5 

5 

45.00 

115.00 

277.5 

9 

43.00 

105.00 

259.5 

5 

47.00 

100.00 

279.5 

10 

46.00 

95.00 

279.5 

6 

46.00 

88.00 

276.5 

10 

44.00 

84.00 

279.5 

6 

43.00 

75.00 

279.5 

5 

46.00 

68.00 

279.5 

9 

34.00 

117.00 

279.5 

5 

33.00 

100.00 

279.5 

6 

33.00 

97.00 

276.3 

2 

32.00 

94.00 

279.5 

4 

35.00 

90.00 

279.5 

5 

32.00 

85.00 

279.5 

5 

38.00 

122.00 

278.5 

6 

41.00 

112.00 

259.5 

1 

43.00 

87.00 

279.5 

1 

42.00 

87.00 

268.5 

1 

36.00 

84.00 

250.0 

1 

41.00 

86.00 

279.5 

9 

38.00 

97.00 

267.5 

6 

35.00 

99.00 

279,5 

5 

38.00 

120.00 

268.0 

1 

35.00 

110.00 

250.0 

1 

45.00 

69.00 

279.5 

1 

39.00 

96.00 

274.5 

1 

40.00 

82.00 

279.5 

1 

40.00 

75.00 

264.5 

1 

40.00 

79.00 

270.0 

1 

35.00 

92.00 

276.5 

1 
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RECOVERY  BASE  DATA 


Latitude 

62.00 

65.00 

65.00 

65.00 

55.00 

52.00 

54.00 

49.00 

64.00 

56.00 

52.00 

60.00 

75.00 

69.00 

65.00 

56.00 

60.00 

52.00 


Longitude 

150.00 

145.00 

157.00 

155.00 

115.00 

107.00 

110.00 

54.00 

68.00 

111.00 

60.00 

125.00 

55.00 

50.00 

20.00 

4.00 

95.00 

-178.00 


PROGRAM  NETWORK 


THIS  PROGRAM  SOLVES  THE  TANKER  TO  SOMBER  TO  RECOVERY  BASE 
MATING  PROBLEM  USINC  A  NETWORK  SOLVER  THAT  MINIMIZES  THE 
TOTAL  TANXER  FLEET  FUEL  CONSUMED.  ONCE  TANKERS  ARB 
ASSICNED  TO  BOMBERS  AND  PRBS,  THE  INDIVIDUAL  REFUEL  INC 
LOCATIONS  ARE  OPTIMIZED  TO  HAZIHIZE  BOMBER  ENTRY  POINT 
FUEL.  BOMBER,  TANXER,  AND  PRB  DATA  ARE  INPUTS  TO  THE 
PROGRAM,  AND  THE  INDIVIDUAL  ASSIGNMENTS  AND  REFUSLINC 
LOCATIONS  ARE  THE  OUTPUTS.  ITERATIONS  CAN  RE  MADE  IN  AN 
ATTEMPT  TO  IMPROVE  THE  ASSIGNMENT  PROCESS.  AND  INDIVIDUAL 
BOMBER  FUELS  AT  THE  ENTRY  POINT  CAN  BE  OUTPUT  AS  WELL. 


COMMON/ FACTOR /PI , RAD 

FILE  DEVICES  USED  IN  THIS  PROCRAN: 


TAPE  1  . BOMBER  INPUT  DATA 

TAPE  2  . TANXER  INPUT  DATA 

TAPE  1  . RECOVERY  BASE  INPUT  DATA 

TAPE  S  . INPUT  FOR  GNET  SUBROUTINE 

TAPE  <  . OUTPUT  FROM  GNET  SUBROUTINE 

TAPE  7  . OUTPUT  FROM  GNET  SUBROUTINE 

TAPE  •  . COST  OUTPUT  FROM  GNET 


INDEI  OF  IMPORTANT  VARIABLES 


NOTRKS . NUMBER  OF  BOMBER  TRACKS 

NOTB . NUMBER  OF  TANXER  BASES 

NOPRB . NUMBER  OF  POST  REFUELING  BASES 

NOTAVL . NUMBER  OF  TANKERS  AVAILABLE 


ARRAYS  SHOULD  BE  DIMENSIONED  AS  FOLLOWS: 


BBLAT ,  BBLOM ,  EPL AT ,  EPLON ,  EPFREQ . NOTRKS «- 1 

niAT.nLON.TKCWT.NOTXRS . NOTB  +  1 

PRBLAT,  PRBLON ,  PRBCAP . NOPRB*  1 

DIET, COURSE, DIFP.UNREF, FUELS, NOTREQ, 

ID  IF,  EPF  ACT,  ZED  IF . NOTRKS 

CPLAT, CPLON , 1ARLAT , EARLCN, FSSTIH, 

ARDI ST , ENDD I B , FUELOF , BONBON , YY , 

FARCVT,  FBI  AT,  PBLON . NOTAVL 

COST1 ,  ARCPGW . NOTB ,  NOTAVL 

COST! ,  EARCWT . NOTAVL ,  NOPRB 

ARDIS, SARDIS . I 
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INTEGER  NOTXRS (41), NOTREQ ( 0 1 > , PRBC AP  <  Z 0 ) 

REAL  IBLAT(ll) , IBLON(ll) .SPLAT (ID , EPLON(Ol) , EPFRSQ(Ol) , 
(TBtAT(41 > , TBLON( 41 ) ,TKGWT(«1 ) , PRBLAT( Zl > , PRBLON< Zl > , 

(DI8T( 0 1 ) , UNREF (ID .FUELS ( II > ,XDIP(I1) .ARDIS(Z) , 

(EARDIS( 3) , COURSE! II ) , CPI ATI  ISO ) >  CPLON( ISO ) , 

(EARL AT( ISO), EARLONI ISO), PESTIK( ISO), COST1 (41,130), 

(C08TZ<  ISO  ,  ZO  )  ,  ARCPGV(41 ,  ISO)  .FUELOPUSO)  , FARCWT!  ISO)  , 

(PBLATt ISO ) ,PBLON( ISO) ,YY(150) ,1XD1P(01) , BONBON! ISO ) , 

(ENDDIS( ISO ) , ARDI ST( ISO ) , EPF ACT ( 11) ,DIPf(tl> , 

(EARCVT! 130,20) 

C 

C  THIS  DATA  STATEMENT  IS  USED  TO  INITIALIZE  ARCP  AND  EAR 
C  LOCATIONS  US1NC  THE  “AVERAGE"  TANKER  CONCEPT. 

C 

DATA  ARDIS , SARDIS/ 1000 . ,3100. , 4000 ., 1030 ., 3243 ., 4120 . / 

PI-3 . 14130Z8S4 
RAD-100. 0/PI 
REWIND  1 
REWIND  Z 
REWIND  3 
REWIND  S 
REWIND  0 
REWIND  7 
C 

READ  THE  INPUT  DATA  FROM  TAPES  1,  Z,  AND  3.  AND  CALCULATE  THE 
NUMBER  OF  BOMBER  TRACKS,  TOTAL  NUMBER  OF  TANKERS  AVAILABLE, 
NUMBER  OF  TANKER  BASES  USED,  AND  THE  TOTAL  NUMBER  OF  POST 
REFUELING  BASES. 

1-1 

10  READ! 1, • ,END-Z0 )BBLAT( I > , BBLON! I ) , EPLATC I ) , EPLON( I ) , EPFREQ( I ) 
1.1*1 
GOTO  10 
SO  NOTRKS-I -1 
1-1 

NOTAVL.O 

30  REAO( S, * , END-40 )TBLAT( I ) , TSLON< I ) , TKGWT( I ) ,NOTKRS( I ) 
MOTAVL-NOTAVL*KOTXRS ( I ) 

1-1*1 
GOTO  30 
40  NOTB-I-1 
1-1 

SO  RBAO! 3, *, END-00 )PRBLAT( I ) , PRBLON< I ) , PRBCAP! I > 

1-1*1 
GOTOSO 

00  NOPRB-I-1 
TAS-444. 

ITER-1 
BIGM-10000 . 

NSOM-O 

FOR  BACN  BOMBER  TRACK,  ASSIGN  THE  APPROPRIATE  NUMBER  OF 
TANXSRS  PROM  THOSE  THAT  ARE  AVAILABLE. 
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SO  <0  1*1 .NOTRKS 

CALL  CIRCLE!B8LAT! I) , BBLON! I) , EPLAT( I ) , EPLONC I ) , X, Y) 

D1ST!I)-Z 

COURSE < I ) «T 

DLOTEF-D I  ST  ( I )  •  1 0  3  . 

TLOTEP- ( DLOTIP/TAS) *10 . 

HOSIC-TLOTEP/ 30. 

TLEFT-TLOTEP-NOSEG*3  0 . 

CVT-470 . 

DO  70  J-l, NOSEC 
ff-io.  n+.  7ii*cw 
IF  <GVT. LE . 340 . )FF-40 . 73* . «ZS»CVT 
FUEL-FF* . 03 
6VT-CVT-PUEL 
70  CONTINUE 

FF-10  Z7+.710*GVT 

IFICVT.LE . 340 . )FF-40 . 73*  .  018*GVT 

GVTEP-GVT-! FF*TLEFT* .001) 

UNREF!I)>GVTEP-Z18 .3 
DIFF( I >-EPFREQ! I ) -UNREF ( I ) 

NOTREQ! I >-l 
FUELS! 1) -03. 

IF(OIFF( I ) . CT. 03 . )THEH 
NOTREQ! I ) -Z 
FUELS! I). 135. 

ENOIF 

IF<D1FF< I) .CT. 133. >TNSN 
NOTREQ! I )»3 
FUELS! I ) -104 . 

ENOIF 

NSUN-NSUH+NOTREQ ! I ) 

00  CONTINUE 
NRIO-O 

I F ( NOTAVL . CE . NSUH ) GOTO  130 
NN I  D-NSUN-NOTAVl 
00  DO  100  J-l, NOTRKS 

XDIF! J) -FUELS! J)-D IFF!  J) 

100  CONTINUE 

FNAE-ZDIF! 1) 

DO  110  J-Z. NOTRKS 

IF!XDIF! J) .  CT .  FHAX ) FMAX-ZDIF ! J) 

110  CONTINUE 

DO  130  I -1. NOTRKS 
IFdDIF!  I  > .  EQ.  FHAZ) THEN 
NRID-MRID-1 
NOTREQ! I ) -NOTREQ! I ) - 1 
IF (NOTREQ! I) .EQ.Z)THEN 
FUELS! I ) -135 . 

ELSE  IF  (NOTREQ!  I)  EQ.  DTHEN 
FUELS! I) -03. 

ELSE 

FUELS! I) -0. 

INDXP 
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GOTO  ISO 
END  IF 

120  CONTINUE 

130  IF(NRID.GT.O)GOTO  IS 
C 

C  DETERMINE  INITIAL  ARCF  AND  EAR  LOCATIONS  FOR  EACH  TRACK 

C  USING  THE  "AVERAGE*  TANKER  DISTANCES.  ALSO  DETERMINE  THE 

C  TIME  REQUIRED  FOR  THE  BOMBER  TO  GET  TO  THE  ARCF.  THIS  TIME 

C  IS  USED  TO  CHECK  FOR  TANXER  FEASIBILITY  AT  EACH  ARCF. 

C 

MN-0 

DO  ISO  I-l.NOTRKS 
DO  140  Jal ,NOTREQ( I) 

NN-NN+1 

CALL  LATLONI BBLAT< I ) , BBLON( I > , ARDIS < J> , COURSE ( I ) , S,T) 

CFLAT(NN)aS 

CFLON(NN)aT 

CALL  LATLONI BBLATt I) , BBLONI I > , EAROI S<  J) , COURSE! I > ,S,T) 
EARLAT(NN) «S 
EARLON(NN) *T 

FESTIM(NN)a( ( ARDIS (J) -104. )*I0 . ) /TAS 
140  CONTINUE 
ISO  CONTINUE 

FRINT* ,'******************************* 

PRINT* 

PRINT*,*  NETWORK* 

PRINT* 

PRINT*, •***••**•**«**•«***•***»•**•*••• 

100  PRINT* <////>* 

PRINT*, '****«*«**•«*«*«*•«•**••••••••*• 

PRINT*,'  FOR  ITERATION  NUMBER  '.ITER 
PRINT*, •«***««*«**••****«*«*•*••******• 

PRINT*  <////!* 

C 

C  DETERMINE  THE  TOTAL  NUMBER  OF  NODES  IN  THE  NETVORK. 

C  AND  WRITE  THIS  AS  THE  FIRST  ENTRY  ON  TAPE  S. 

C 

HaNOTB+NOTAVL *  2  ♦NOPRB* 1 
WRITE (3,2S0)H 

DETERMINE  ARC  COSTS  FROM  EACH  TANKER  BASE  TO  EACH  ARCF, 

AND  WRITE  THIS  ON  TAPE  S  ALONG  WITH  THE  UPPER  AND  LOVER 
BOUNDS  FOR  EACH  ARC. 

C 

DO  200  I  a 1 , NOTE 
DO  100  Jal  ,JOi 

CALL  CtRCLEITBLATd) .TBLOMI I),CPLAT< J> ,CPL0N( J) ,1, Y> 

Dismcai 

TIMER*! (DISTNC-104 . ) *00 . > /TAS 
NSEC-TIMER/30. 

TL  EFTaTI MER-NS EG  *  1 0 . 

CVaTXCVTI I ) 

DO  170  Xal.NSEG 
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FFa . 83I*GW»2 . II 
IF(CW. LE . 180 . )FF> . 813*CV«10 . 8 
FUCL>FF* .01 
GVbGV-FUIL 
170  CONTINUt 

FF-.S30*GW*2.83 

IF  (GW.  H.  180.  )FFa.  8  1  J*GW+10  .  8 

CWaCW- < ( FF'TLEFT) *001) 

C 

C  DETERMINE  IF  HOLDING  18  REQUIRED 
C 

FESTIH<J>-FESTIM<J)+.001 
IF (TIMER . CT. FSSTIM< J  > )THEM 
C08TKI,  J)-IIGM 

ELSE 

HTIHEaFESTlMU) -TIMER 
M8SG-HTIME/30. 

TLEFTaHTIME-MSEG*  30 . 

DO  180  Xal.NSEC 
FFa . 0023*GW- 18 . 18 
CVaGW-FF* . 03 
180  CONTINUE 

FFa . SO  SS*CW-18 . 88 
CWbGW-FF*TLEFT/ 1000 . 

ARCPCWI I , J) -CV 
C08TKI,  J)-110.  1-GW 
END  IF 

ICOST-COSTKI.  J>*10000 
NODES  aMOTI ♦ J 

WRITE ( 3 . 270 ) I .NODES . I  COST , NOTKRS ( I ) , 0 
100  CONTINUE 
200  CONTINUE 
C 

C  WRITE  REFUELING  ARC  DATA  ONTO  TAPE  3. 

C 

ICOSTaO 
DO  210  lal.NN 
NODEaNODES* I 
JaNOTB+I 

WRITE ( 3, 270 >J, NODE, I  COST, 1,1 
210  CONTINUE 

DETERMINE  ARC  COSTS  FROM  EACH  EAR  NODE  TO  EACH  PRB  NODE. 
THIS  ARC  COST  INCLUDES  A  PROVISION  FOR  NOT  CLIMBING  TO 
C  OPTIMUM  ALTITUDE  IF  THE  PRB  IS  NEARBY.  WRITE  THE  COSTS 

C  AND  UPPER  AND  LOVER  BOUNDS  ONTO  TAPE  3. 

C 

FRBCWTal IS . 1 
DO  280  Ial.MI 
DO  280  Jal.NOPRB 

CALL  CIRCLE  (EARL  ATI  I ) ,  SARLONI I  > ,  PRBLATU) ,  FRBLONt  J) ,  E,  Y ) 

CW-PRBCWT 

IF ( X . LE . 180. >TMEN 
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TLKFTaI/TAS*00. 

GVaCW+ . 197*TLBFT 

I  LSI 

DISTXaX-73 
TXHEaOISTX/TASMO . 

NSECaTlMK/30. 

TLEFT.TIHE-NSXC*  30 . 

00  220  Kal.NSEG 

CVa  <0V+ .3241) / . 9731 
220  CONTINUE 

GVa6V*<< . 01 3* ( CV+ 1 . 3 ) +10 . 8  > *TLEFT/ 1000 . )*2. 

wait 

EARCWTd,  J)-CW 

COST2! I ,J)aEAHCWT(I , JJ-PRBCVT 
I COSTaCOSTZ ( I , J ' *10000 
INODESaNODE+J 
I IaNODES+1 

WR ITE ! 9 , 270  > 1 1 , INODES , I  COST, 1 , 0 
230  CONTINUE 
240  CONTINUE 
C 

C  WRITE  ONTO  TAPE  S  THE  ARC  DATA  FROM  EACH  OF  THE  FRBS 
C  TO  THE  SECONDARY  SINK  NODE. 

C 

ICOSTaO 

DO  290  Ial.NOFRE 
JJaNODE+I 
I IaINODES+1 

WRITE (3,270)JJ,II,I COST , PRBCAP! I ) , 0 
290  CONTINUE 
C 

C  WRITE  ONTO  TAPE  9  THE  ARC  DATA  FROH  THE  SECONDARY 
SI NX  TO  THE  SUPER  SINX. 

JJ-IW2 

WRITE (9, 270) 1 1, JJ, ICOST,NN,0 

WRITE  ONTO  TAPE  9  THE  ARC  DATA  FRON  THE  SOURCE 
NODE  TO  EACH  OF  THE  TANKER  BASE  NOOES. 

DO  200  lal.NOTB 

JJalUl 

WRITE (9, 270) JJ,  I,  ICOST,NOTKRSd ) ,  0 
200  CONTINUE 

270  FORMAT! 11,211,22,3110) 

200  FORMAT! 19) 

200  FORNAT!9S, 219, 9110) 

300  FORMAT! IE, F12.0) 

REWIND  9 

CALI  THE  CNST  SUBROUTINE  TO  SOLVE  FOR  THE  MINIMUM  COST 
FLOW  THROU6H  THE  NETWORK  THAT  HAS  BEEN  DESCRIBED  BY 
TAPE  9.  THE  ENTIRE  CHET  OUTPUT  IS  STORED  ON  TAPE  0, 
AND  THE  ARC  FLOWS  ONLY  ARE  STORED  ONTO  TAPE  7. 
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CALL  GNETSS 
REWIND  I 

READ!  •>  3O0)RCOST 

PRINT* , *  INITAL  COST  IS  * .RC08T 
PRINT’ (///) ’ 

REWIND  7 
C 

C  READ  TAPE  7  AND  DETERMINE  THE  PR8  MATINGS  FROM 

C  EACH  OP  THE  EAR  POINTS.  ALSO  DETERMINE  THE  TANKER 

C  CROSS  WEIGHTS  AT  EACH  OF  THE  EAR  POINTS. 

C 

NTIHESaNOTB*NOTAVL+NOTAVL 
DO  310  IAal , NTIME3 
READ! 7 ,210) IT.NH. ISL , IX , I  CP , IC , IRC 
310  CONTINUE 

DO  330  I  Sal , NOPRB 
DO  320  ISal.NOTAVL 
I BBaNOTAVL+NOTB+ I S 
IXCa2*NOTAVL*NOTB+ II 
READ(7,2tO>IT,MH,IBL,IX,ICP,IC,IRC 
IF<  IX  .  GT.  0  .  AND.  IT.  EQ .  IBB  .  AND  .  NH .  EQ.  IXOTHEN 
FARGWTI IS) aE ARCWT  (IS, IB) 

PBLAT( IS) aPRRLATC I B) 

FBLON ( 18 ) aPRBLONC II) 

END  IF 

320  CONTINUE 
330  CONTINUE 
REWIND  7 
C 

C  READ  TAPE  7  ACAIN  TO  DETERMINE  THE  TANKER  TO  BOMBER 
C  ASSIGNMENTS.  ALSO  DETERMINE  THE  FUEL  CONSUMED  ON 
C  THE  AR  TRACK  BT  BOTH  THE  BONIER  AND  TANKER.  DETERMINE 

C  TANKER  OFFLOAD  CAFAIILITT  AND  BOMBER  (BILOAD  CAPABILITY 

C  FOR  THE  MATINGS  GIVEN. 

C 

DO  300  Jal.NN 
DO  3B0  Ial.NOTB 

READ( 7, 200 > IT,NH, IIL , IX , ICP, IC , IRC 
JJ-NOTB+J 

IF< IX . GT. 0 . AND . I . EO. IT. AND. JJ . EQ . NH)THEN 

CALL  CIRCLE(CPLAT( J) , CPLON( J ) , EARLATI J) , SARLON( J ) ,X , Y) 
DISTNC.X 

TIMEaDISTNC/400 . *00 . 

FPa. 524* ( ( ARCPGWI I , J) +FARGWT( J) ) /l) +07 . 71 
FUXl-FF*TIHE/ 1000 

FUELOF ( J ) a ARCPCWI I , J ) -F ARGWT( J > -FUEL 
MMaO 

DO  370  Kal.NOTRKS 
DO  300  Lal.NOTREQ(X) 

MM.HM ♦  ! 

1F< J . EQ . MM)  THEN 
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CALL  CIRCLE! BBLAT( X) , BBLON! K) , CPLAT!J ) , CP LON < J) , X, Y > 

DISTY. I 

ARDIST!HM).X 

TIMER. ( (DUTY- 10 3  .  >/TAS)«IO. 

NSEC. TIMER/  30. 

Tt EFT. TIMER -NSEC* 30 
CV.470 . 

00  340  H-l.NSEC 

CV.CV- ((10.27+  7 18*GV) * . 03) 

340  CONTINUE 

BCWTCP.CV- ((10. 27*. 7 1 8*CV) *TLEFT/ 1000 .  > 

IFtL.CT. DTHEX 

CALL  CIRCLE! EARLAT( J- 1 ) , EARLON! J-l ) , CFLAT( J) , 
t  CPLON(J) ,Z , Y) 

OISST.X 

TIMES- DISST/TAS*80. 

NSEC. TIMES/ 30. 

TLEFT.TIMES-NSEC*30 . 

CW.488 . 

OO  390  IT.l , NSEC 

CV.CV- ( < 10 . Z7«  . 710*CV) *  .  03) 

350  CONTINUE 

BGVTCP.CV-! <10 . tl+ . 718*GV> *TLEFT/ 1000 . ) 

END  IF 

FUEL-TINE* <80. 93+ <BGVTCP*488. >  /  2* . 78  37) / 1000 . 

BONBON  <  J ) .4  8 1 . - <  BGVTCP-FUEL) 

END  IF 

380  CONTINUE 
370  CONTINUE 
END  IF 

380  CONTINUE 
380  CONTINUE 
C 

C  ARCPS  AND  EARS  ARE  NOV  ADJUSTED  A  CERTAIN  DISTANCE 
C  IN  AN  ATTEMPT  TO  EQUATE  BOMBER  ONLOAD  CAPABILITY  WITH 
C  TANKER  OFFLOAD  CAPABILITY.  THE  DISTANCE  THE  ARCP  IS 
C  MOVED  DEPENDS  ON  HOW  MUCH  THESE  TWO  NUMBERS  DIFFER 

C  INITIALLY.  NEW  BOHBER  TIMES  TO  THE  ARCFS  ARE  ALSO 

CALCULATED. 

400  J.O 

DO  420  I -1 , NOTRKS 
DO  410  X.l.NOTREQd) 

J.J»1 

X. FUELOF <  J ) -BONBON! J ) 

IF !ABS!X)  CT. . 4ITHKN 
IFIX.CT. 0)THEN 

IFIX.LT. . 3)THEN 
ARDIST! J).ARDIST( J)+l . 

ELSKIFIX.LT. 1 . )THEN 
ARDIST! J).ARDIST( J)+10  . 

ELSEIFIX.LT. 2. ) THEM 
ARDIST! J).ARDIST( J)*13 . 
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ELSEIF <Z. IT. 3 . I THEN 
ARDISTt J) .ARDISTt J) +30  . 

ELSEIF (Z . IT. 4. )THEN 
ARDISTt JI. ARDISTt JI+40 . 

EL3EI F  <X . IT. 7 . ITHEN 
ARDI3T( J I. ARDISTt JI  +  83 . 

ELSE  I F ( Z . LT . 1 0 . 1 THEM 
ARDISTt  J I .ARDI 3T( J) +100 . 

ELSE  I F ( X . LT . 1 S . ) THEN 
ARDISTt J > «ARDI STt J I +200 . 

EL3EIFd.LT.  30. ITHEN 
ARDISTt J) -ARDISTt JI+300. 

ELSE  IF ( Z . GT. 30. )THEH 
ARDISTt J) •ARDISTt J) +400 . 

END  IF 

ELSEIFtE . LT. 0 . ITHEN 

IF  t ABStZ) .LT.  . 3) THEN 
ARDISTt JI .ARDISTt J) -4. 

ELSEIF (ABStZ )  IT. 1. ITHEN 
ARDISTt J) .ARDISTt J) -8. 

ELSEIF (ABStZ) . LT . 2. ITHEN 
ARDISTt  J) .ARDISTt  J) -14 . 

ELSEIF (ABStZ)  LT. 3. ITHEN 
ARDISTt J). ARDISTt JI-28 . 

ELSEIF (ABStZ) .LT. 4. ITHEN 
ARDISTt J I .ARDISTt J) -38. 

ELSEIF (ABS(Z) .LT. 7. ITHEN 
ARDISTt JI.ARDl STt J I -82. 

ElSEIF(ABSd)  LT.  10.  ITHEN 
ARDI STt J I .ARDI STtJ >-95. 

ELSEIF  <  ABStZ) . LT. IS . ITHEN 
AFDIST(J). ARDI STtJ) -185. 

ELSEIF (ABS(Z) .LT. 30 . ITHEN 
ARDISTt J). ARDISTt JI -280. 

ELSEIFt ABStZ) . CT. 30 . ITHEN 
ARDISTt J I- ARDI STtJ) -380. 

ENDir 

END  IF 
END  IF 

ENDOIStJ). ARDISTt J) +130 . 

IFtK. GT. 1 1 ENDDIS(J). ARDISTt J) +  83. 
irtK.CT. ZIENDDISt J) .ARDISTt JI+80 . 

CALL  LATLONt  BBLATt 1 1 , BBLONt 1 1 .ARDISTt J) , COURSE ( I 

CPLATtJI.S 

CPLONtJI.T 

CALL  LATLONt BBLATt I) .BBLONt I) .ENDDISt J) , COURSE ( I 
EARLATt J 1.8 
EARLONt JI.T 

FESTINt J)«t  t ARDISTt J) -184. 1 *80. I /TA8 
410  CONTINUE 
420  CONTINUE 


I.S.T) 

I.S.T) 
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C  THE  REMAINDER  OF  THE  PROCRAH  RECOMPUTES  THE  TANKER 
C  OFFLOAD  CAPABILITY  AND  BOMBER  ONLOAD  CAPABILITY  AND 
C  REITERATES  THIS  PROCEDURE  UNTIL  THE  TWO  NUMBERS  ARE 

C  WITHIN  400  POUNDS  OF  EACH  OTHER.  BOMBER  ENTRY  POINT 

C  FUEL  IS  COMPUTED  FOR  EACH  BOMBER,  AS  WELL  AS  THE  TOTAL 

C  ENTRY  POINT  FUEL  FOR  THE  BOMBER  FLEET.  ALSO,  THE  DIFFERENCES 

C  BETWEEN  BOMBER  FUEL  REQUIRED  AND  ACTUAL  BOMBER  FUEL  IS 

C  CALCULATED  FOR  EACH  BOMBER. 

C 

00  440  Jal , NOTAVL 

CALL  CIRCLE(PBLAT( 4) , PBL0N< J> , EARLATt J> ,EARLON< J> ,X, Y) 
GWsPRBGWT 
IPtE.LE. 140. )THEN 
TLEFTaZ / TAS*  0  0 . 

GW.CW*.157*TLEFT 

ELSE 

DI STXaX-73 . 

TIMEaDISTX/TASMO . 

NSECaTIME/ 30 . 

TLEFT-TIME-NSEG* 30 . 

DO  430  KXal.NSEG 

CWa(CW+ . 3241 ) / . 07SI 
430  CONTINUE 

GW.GW><<  013*<GW«1.S>«10.0)*TLEFT/1000. >♦!. 

END  IF 

F  ARGWT (J) aCW 
440  CONTINUE 
REWIND  7 

DO  310  Ial, NOTAVL 
DO  3X0  J  a I , NOTB 

RE AD (7, 110) IT.NH, IBL, IX, ICP, IC, IRC 
JJaNOTB+I 

I F  < IX. GT . 0 . AND . J . EO . IT. AND . JJ . EQ . NH)THEN 

CALL  CIRCLE (TBLAT(J) , TBLON ( J ) , CPLAT ( I ) , CPLONI I ) ,S ,T) 
DISTT.S 

TIMER- (DISTT- 104 . MOO. /TAS 
NSEGaTIMER/  30 . 

TL EFT -TIMER -NS EG*  30 . 

GW-TXGWT(J) 

DO  430  Lai ,NSEG 

FFa . 030*GW+Z . 03 
IF (GW. IE. 100. >FFa.013«CW*10.0 
GWaCW-FF* . 01 
430  CONTINUE 

FFa . S30*GW+t . 03 

IF(CW.LE. 100. IFF-. S1J*CV^10.I 

CW-CW-FF*TLEFT/ 1000 . 

HTIME-FESTIH(  I ) -TIMER 
NSEG-HTIME/30. 

TLEFT-HTIME-NSEC* 10 . 

DO  400  Hal, NSEC 

FFa.00X3oGW-lB.00 
GWaCW-FF* .  03 
CONTINUE 


400 
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FFa.9023«GV-18.89 
CV.GV-FF*TLEFT/1000. 

ARCPCVC J , I )«GW 

CALL  CIRCLE ( CPLATC 1 ) , CPLON( 1 > , EARLATC I > , EARLON( I > , I , Y> 
DISTTaX 

TIMEaDISTT/400 . »80 . 

FFa . 3  24*  <  ( ARCPGVC  J , 1 > ♦FARGVT< I  >  > /!) 487 . 71 
FUELaFF*TINE/ 1000 . 

FUELOF ( I ) aARCFCVf  J,I>- FARGWTI I > -FUEL 

rauo 

DO  300  X-l,NOTRKS 
DO  490  L.l ,NOTHEQ(X> 
fw-m+i 

!F( I . EQ . MM) THEM 

CALL  CIRCLE ( BBLATC  X)  , BBLON(K) . CPLATt t ) . CPLONI I ) . X , Y> 
DUTY. I 

TIHER.((DISTY-105 . >/TAS>*80. 

NSECaTIMER / 30 . 

TLEFTaTIHER-NSEC*30 
GV.479. 

DO  470  Hal ,N3EG 
GW-CW-< (10 . 27* . 718*CV>* .03) 

470  CONTINUX 

BGVTCP.GV- ( ( 10.274. 718»CW) *TLEFT/ 1000 . > 
IF(L.CT.1>THEN 

CALL  CIRCLE (EARL ATI  1-1  >  # EARLONI 1-1)  , CPLAT( I > , 

A  CPLON< I ) , X , Y ) 

DI88T.X 

TIME8-DI88T/TA8*«0. 

N8EG. TIMES/ 30. 

TLEFTaTINBS-NSEG*30 . 

CVa48t. 

DO  490  ITal ,NSEG 

CVaGV- ((10.274.7 1 l*GV) *.03) 

490  CONTINUX 

BCVTCpaGV- ( { 1 0 . 274. 718«CW)*TLEFT/ 1008. ) 

END  IF 

FUEtaTINE* ( 90 . 334 (ICVTCP4489 . ) / 2. ». 7837 ) / 1800 . 
BONBON ( I ) .499 . - ( BGVTCP-FUEl) 

END  IF 

490  CONTINUE 

300  CONTINUE 

ENDIF 

310  CONTINUX 
3X0  CONTINUE 

DO  330  Ial.NOTAVt 

FT ( I) aAB8(FUBL0F( I > -BONBON ( I >  > 

330  CONTINUE 
TTTaTY(l) 

DO  340  JaX.NOTAVl 

IF(TT ( J) . GT . TTY)YTT«TT( J> 

340  CONTINUE 

IFIYYY.CT. . 4 > GOTO  400 


fsum-o. 

JaO 

DO  SiO  1-1 , NOTRKS 
J-J+NOTREQC I ) 

IF (NOTREQI I )  EO.OITHEN 
EPFACTI X  laUNRSF ( I > 

USE 

CALL  CIRCLE! EPLATt I > , EPLON< I > , EARLATC J) , EARLON! J) , I , Y) 
DISTT-X 

TIMERaDt  STT/TASMO . 

NSEC -TIMER/ 30. 

TL  EFT-TI MER -NS  EC  *  3  0 . 

CV-48B . 

DO  SSO  Xal.NSEG 

IF (GW.CT. 340. >FFa 10.37*. 71MGW 
I F ( CV .  LE  .  340.  >FFa40.73+  .m*CV 
GWaCW-FF* . 03 
330  CONTINUE 

IF(CV.CT.340.  )FFalO.Z7*.7H*CW 
I F  <  CV. LE . 140 . )FF-4# . 73* . IZ8*CV 
CVaCW-FF*TLEFT/ 1000 . 

EFFACT!  DaCW-tlO .  I 
END  IF 

FSUHaFSUM+EPFACT( I > 

EXDIF ( I ) -EPFACT! I ) -EPFHEQ  < I ) 

300  CONTINUE 
REWIND  3 
REWIND  • 

REWIND  7 
REWIND  0 
ITXRaITER+1 

IF ( ITER . LE . 10) GOTO  110 
S' 


n  n  n  n  n  n 


PROGRAM  GREED? 


C 

C  THIS  PROGRAM  SOLVES  THE  TANKER  TO  BOMBER  TO  RECOVERY  BASK 
C  MATING  PROBLEM  USING  A  "CREED?"  TYPE  ALGORITHM.  ASSIGNMENTS 

C  ARE  HADE  BASED  ON  COST  DIFFERENCES  BETWEEN  THE  TWO  BEST 

C  POSSIBLE  ASSIGNMENTS  FOR  A  GIVEN  REFUELING  LOCATION,  AND  THE 
C  REFUELING  WITH  THE  BIGGEST  DIFFERENCE  IS  ASSIGNED  FIRST. 

C  AFTER  THE  ASSIGNMENTS  ARE  MADE,  INDIVIDUAL  REFUELING  LOCATIONS 
C  ARE  OPTIMIZED  TO  MAXIMIZE  BOMBER  ENTRY  POINT  FUEL.  BOMBER. 

C  TANKER,  AND  PRB  DATA  ARE  THE  INPUTS  TO  THE  PROGRAM.  AND  THE 
C  INDIVIDUAL  ASSIGNMENTS  AND  REFUELINC  LOCATIONS  ARE  THE  OUTPUTS. 

C  ITERATIONS  CAN  BE  MADE  IN  AN  ATTEMPT  TO  IMPROVE  THE  ASSIGNMENT 

C  PROCESS,  AND  INDIVIDUAL  BOMBER  ENTRY  POINT  FUELS  CAN  BE 
C  OUTPUT  AS  WELL. 

C 

C 

COMMON/ FACTOR/ PI , RAD 
C 

C  FILE  DEVICES  USED  IN  THIS  PROGRAM: 


C 

C  TAPE  l  . BOMBER  INPUT  DATA 

C  TAPE  2  . TANKER  INPUT  DATA 

C  TAPE  3  . RECOVERY  BASE  INPUT  DATA 

C  TAPE  7  . TANKER /BOMBER /PRB  ASSIGNMENTS 

C 

C 


C  INDEX  OF  IMPORTANT  VARIABLES 


C 

C 

C  NOTRKS . NUMBER  OF  BOMBER  TRACKS 

C  NOTE . NUMBER  OF  TANKER  BASES 

C  MOPRB . NUMBER  OF  POST  REFUELINC  BASES 

C  NOTAVL . NUMBER  OF  TANKERS  AVAILABLE 

C 

C 


C  ARRAYS  SHOULD  BE  OIMINSIONED  AS  FOLLOWS: 


C 

C  BBLAT,  BBLON .  EPl AT ,  EPLON ,  EPFREQ . NOTRKS*  1 

TBLAT ,  TBLON ,  TKCWT ,  NOTXRB ,  NOTXS . NOTB+1 

PRBLAT,  PRBLON,  PRBCAP ,  PRBCP . NOPRB+1 

DIST, COURSE, DIFF.UNREF, FUELS, NOTRSQ, 

XDIF.EPFACT,  XXDIF . NOTRKS 

CPLAT , CPLON.EARLAT, EARLON , FISTIH, 

ARD I ST , ENDD I S , FUE  LOF , BONBON , YY , 

C  BEST1,  BESTS . NOTAVL 

C  COST1 ,  ARCPGW . NOTE ,  NOTAVL 

C  COSTS,  EARCWT . NOTAVL , MOPRB 

C  COSTS . NOTAVL ,  NOTE ,  NOPRB 

C  ARDIS ,  EARDIS . I 
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INTEGER  N0TKR8( 33 ) ,NOTREO( SI ) , PRBCAPI 20 ) 

(, NOTES (33) , PHBCF(  20) 

REAL  BSLAT( II ) , BBLON( II ) , EFLAT( (1 > , EPLON( II >  > EPFREQ( (1  > , 
(TBLAT( 31) , TBLON( 33 ) ,TXCVT( 13 ) , PRBLAT( 20 ) ,  PRBLON( 20  > , 

(DIST( (1 > , UNREF ( (1 > .FUELS (ID ,XO!F(ll) , AR0I5C I) . 

(SARD IS ( 3 ) , COURSE ( II >  >  CPLAT( 133), C?LON< 133 ) . 

(EARL ATI  135 ) ,  EARLONI 13S ) , FESTIHI 113  >  , COST1 (41,135), 

(COST2 (135.20) , ARCPCV( 33,133) , FUELOF ( 113 ) , 

(YY(133>  ,SEDIF((1)  ,BOMBON(13S> , 

(EMDDIS( 133  > , ARDIST( 135) , EPFACT( (1) ,DXFF(133>, 

(S ARCVT( 133,20), COST3 (133,33,20), BEST1 (133), BEST2 ( 133 ) 

C 

C  THIS  DATA  STATEMENT  IS  USED  TO  INITIALIZE  ARCP  AND  EAR 
C  LOCATIONS  USINC  THE  "AVERAGE"  TANKER  CONCEPT. 

C 

DATA  ARDIS, EARDIS/ 1(00 . ,31(0. ,4000. ,1130. ,3243. ,4120. / 

PUS  .  141312(34 
RAD* 1(0 .0/PI 
REVIND  1 
REWIND  2 
REVIND  3 
REVIND  7 
C 

C  READ  THE  INPUT  DATA  FROM  TAPES  1,  2,  AND  3,  AND  CALCULATE  THE 

C  NUMBER  OF  BOMBER  TRACKS,  TOTAL  NUMBER  OF  TANKERS  AVAILABLE, 

C  NUMBER  OF  TANKER  BASES  USED,  AND  THE  TOTAL  NUMBER  OF  POST 
C  REFUELING  BASES. 

U1 

10  READ( 1 . • , END-20 ) BBLAT( I ) , BBLON( I > . EFLAT( I ) , EFLON( I ) . EPFREQI I ) 
Mil 
GOTO  10 
2(  NOTRKB-I-1 
1-1 

NOTAVL-O 

(0  REAO( 2. • .END-40 )TBLAT( I ) , TBLON ( I ) , TKCVT ( I ) , NOTKRS( I > 

NOTAVL -NOTAVL ♦NOTKRS ( I ) 

1-1*1 
GOTO  30 
40  NOTB-1-1 
1-1 

SO  REAO( 3, *, END-40 ) PRBLAT( I ) , FRBLON( I > , FRBCAP ( I ) 

1-1*1 
COTO30 
(0  NOPRB-I-1 
TAB-444. 

1TEH-1 
BICH-((S((. 

NSUM-0 


FOR  EACH  BOMBER  TRACX,  ASSIGN  THE  APPROPRIATE  NUMBER  OF 
TANXSRS  FROM  THOSE  THAT  ARB  AVAILABLE. 


DO  10  I«1 , NOTRKS 

CALI  CIRCLE! BBLAT( 1 ) , BBLON< 1 ) . BPLAT! I ) , EPLON! 1 ) ,I.Y) 
DIRT! I  >aX 
COURSE! I) at 
DLOTEPaDIST! 1 )-105. 

TLOTEPalDLOTEP/TAS) *00 . 

NOSEG-TLOTEP/ 30. 

TLEFTaTLOTEP-NQSCG*SO . 

CWT»470 . 

00  70  Jal, NOSEC 
FFalO . 11* .  71»*CWT 
1F(GUT.LE. 340 . )Ff«40 .73+ . «I»»CWT 
FUBLaFF* .03 
CVTaCVT-FUEt 
70  CONTIMVE 

FFaI0.t7*.710«CVT 

IF (CVT. LE . 340 . >FF«40 .73* . Oil'CVT 

CWTEP«CVT-(FF*TtEFT* . 00 1 > 

UNRIP ( I ) »CWTBP -310  3 
D I FF ( I > aEPFREQ! 1 ) -UNRBF! I ) 

NOTREQ!  Dal 
FUELS! I ) -03 . 

IF (DI FF  < I ) . GT. 0  3 . 1THEN 
NOTREQ ( I )  aZ 
FUELS! I)al 33. 

END  IF 

IFIDIFFd)  -CT.135.  )THBN 
NOTREQ! I ) al 
FUELS! I >  *104 ■ 

END  IF 

NSUHaNSUN*NQTREQ ! I ) 

•0  CONTINUE 
NR 1 Da  0 

IF(NOTAVL.GE  NSWOCOTO  130 
NR I DaMSUM-MOT AVI 
•0  DO  100  Jal.NOTRES 

XDIF! JJaFUELS! JJ-DIFF ! J> 

100  CONTINUE 

FNAXaXDIP! 1) 

DO  110  JaX(NOTRKS 

IFIEDIF! J) . 6T . FNAX ) FMAEaEDI F ! J) 

110  CONTINUE 

DO  1E0  lal  , NOTRKI 
I F  <  XD I F  <  I )  .  EQ . FNAX  >  THEN 
NXID-NRID-1 
NOTREQ ( I ) aNOTRBQ ! I ) - 1 
IF!NOTRXQ!I) . SO. X)THEN 
FUELS! I >alJ3 . 

ELSE IF (NOTREQ! I) . EQ. 1 JTHEN 
FUELS!  DaOI. 

ELSE 

FUELS! I ) aO . 


i 


END  IF 
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GOTO  130 
END  IF 

120  CONTINUE 

130  IF (NRID. GT. 0>  GOTO  00 
C 

C  DETERMINE  INITIAL  ARCP  AND  EAR  LOCATIONS  FOR  EACH  TRACK 

C  USING  THE  “AVERAGE"  TANKER  DISTANCES.  ALSO  DETERMINE  THE 

C  TINE  REQUIRED  FOR  THE  BOMBER  TO  GET  TO  THE  ARCP.  THIS  TIME 

C  IS  USED  TO  CHECX  FOR  TANXER  FEASIBILITY  AT  EACH  ARCP. 

C 

NNaO 

DO  ISO  Ial , NOTRKS 
DO  140  Ja 1 , NOTREQ ( 1 > 

NNaNN+l 

CALL  LATLONI BBLAT< I ) , BBLON( I > , ARDISI J> , COURSE! I ) ,S ,T> 

CPLAT(NN) aS 

CPLON(NN)aT 

CALL  LATLONI BBLATI I ) ,  BBLON( I > , EARDIS ( J> , COURSE! I )  ,  5 , T) 

EARLAT(NN)aS 

EARLON(NN)aT 

FESTIM(NN) a( ( ARDIS ( J) -114 . >  *10 . ) /TAB 
140  CONTINUE 
ISO  CONTINUE 

PRINT*, ••*«••«**••**•*•**•****•••***••* 

PRINT* 

PRINT*, •  GREEDY* 

PRINT* 

PRINT*, •«•••**••**«•****•*•*•*****•**•• 

140  PRINT* <////)* 

PRINT*, * «****•*«•••*•••*•••••*••••••*• • 

PRINT*,’  FOR  ITERATION  NUMBER  *,ITEB 
PRINT*, ******************************** 

DO  101  Ial.NOPRB 
PRBCP ( 1 ) aPRSCAP ( I ) 

111  CONTINUE 

DO  IfE  Ial, NOTE 
NOTKS ( I ) aNOTXBS ( I ) 

102  CONTINUE 
C 

C  DETERMINE  COSTS  FROM  EACH  TANXER  BASE  TO  EACH  ARCP 
C 

DO  200  lal.NOTB 
DO  100  Jal.NOTAVl 

CALL  CIRCLE(TBLAT( I ) ,TBLON( I ) , CPLAT( J) ,CPLON(J> ,E,Y) 
DISTNCaX 

TIMERal (OI STNC-114 . ) *10 . ) /TAS 
NSEG-TIMER/10. 

TLEFT.TIMER-NSEC«30 . 

CVaTXCVT ( I ) 

DO  170  Xal.NSEG 

FPa . 0SI*CV«2 . 03 

IFICV.LE. 100. »rF-.«ll*CV*10.l 

FUELaFF* . 01 

GVaCV-FUEL 
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170  CONTINUE 

TTm . S3I*CV*3 . 03 

IF (CU. LB. 100. >FF..313»GV*10.0 

GVaGV- < ( FF*TLEFT) *.001) 

C 

C 

C  DETERMINE  IF  HOLDINC  IS  REQUIRED 

C 

C 

FESTIMI J)»FESTIN< J>* .001 
IFITIMER.CT. FESTIM(J) >THEN 
C0ST1 ( I , J) aBICH 

ELSE 

HTIHE-FESTIH<J> -TIMER 
NSECaHTINE/ 30. 

TLEFT«KTIME~NSEG*30 . 

DO  ISO  Kal ,N8EC 
FFa . 002S*CV-1S . 00 
CWaCW-FF* . 03 
ISO  CONTINUE 

FF-.00Z3*GV-1S . 01 
CW-CW-FF  *TLEFT/ 1000 . 

ARCPGV! I i J) aCV 
COST1 < I <  J)al 10 . 1-CW 
ENDXF 

100  CONTINUE 
300  CONTINUE 
C 
C 

C  DETERMINE  COSTS  FROM  EACH  EAR  POINT  TO  EACH  PRB 
THIS  COST  INCLUDES  A  PROVISION  FOR  NOT  CLIMBING 
TO  OPTIMUM  ALTITUDE  IF  THE  PRB  IS  NEARBY 
C 
C 

PRBOVTa 113.1 
DO  300  lal  ,  NOT  A  VI 
DO  330  Jal.NOPRB 

CALL  CIRCLE! EARL ATI  I ) , EARLONI I > , PRBLATI J> , PRBLON! J) , E, Y) 
CVaPRBGVT 
IP(E. IE. 140 . (THEN 
TLSFTaS / TAB*  0  0 . 

GVaCVa . 137*TLEFT 

ELSE 

DISTEaS-73 . 

TIHEaDIST2/TAS*00 . 

MSEC-TIME/ 30. 

TLEFTaTIME-NSEC* 30 . 

DO  330  KX>1,NSEC 

CV-(CV«. 3341)/. 0730 
330  CONTINUE 

6Va6V*(< .013*<CV*1.S>*10.0>*TIEFT/1000. )*3. 

END  IF 

EARGVT! I , J)aCV 

COST3! I , J>  aEARGVTI I , J) -PRBCWT 
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130  CONTINUE 
140  CONTINUE 
C 
C 

C  DETERMINE  TOTAL  COST  rOR  EACH  TANKER  BASE, 

C  REFUELING  LOCATION,  PRB  COMBINATION.  THIS 

C  IS  DONE  ST  AODINC  TOGETHER  THE  TWO  COSTS 
C  PREVIOUSLY  CALCULATES. 

C 

C 

DO  ISO  lal , NOTAVL 
DO  130  Jal.NOTB 
DO  143  Xal.NOPRB 

COST! ( I  •  J ,  X ) aCOSTl ( J , I > ♦COST! ( I , X) 

IF (COST J< I , J,K) . CT. RICH) COST1 < I , J , K> -BICM 
143  CONTINUE 
130  CONTINUE 
100  CONTINUE 
NN-NOTAVL 
C 

c 

C  FIND  THE  TWO  BEST  COSTS  FOR  EACH  REFUEL  INC, 

C  AND  FROM  THESE  FIND  THE  REFUELING  WITH  THE 

C  LARGEST  DIFFERENCE  AND  BEST  COST.  THIS 
C  REFUELING  VILL  BE  ASSIGNED  FIRST.  WRITE 
C  THE  REFUELING  ASSIGNMENTS  ONTO  TAPE  7. 

C 

C 

170  BIGDIF-0. 

DO  100  lal .NOTAVL 
BRSTKDaBlGH 
BESTS  ( I )  all  CM 
DO  100  Jal.NOTB 
DO  100  Xal.NOPRB 
IT  icon*  (l.O.K) .It.  BEST!  <  I  >  JTHEN 
BESTS ( I ) -COSTS ( I , J , K) 

KFACEaJ 
END  IF 

100  CONTINUE 
100  CONTINUE 

DO  SOI  Jal.NOTB 
DO  SOS  Xal.NOPRB 

IF (COSTS! I . J , K) . GE . BIST1 ( I > . AND . COSTS ( I . J. X)  LE. BESTS (I)  AND. 
IKFACE .NX.  J) BESTS ( I ) -COSTS < I , J , X) 

SOt  CONTINUE 
SOI  CONTINUE 

DIPF( I >  aBESTS ( I > -BEST1 ( I ) 

1F(DIFF(I)  GT.BICDIFITHEN 
BICOI FaDIFF < I ) 

IFACE-I 
END  IF 

SOO  CONTINUE 

DO  SOO  Jal.NOTB 
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00  330  Kal , NOPRB 

IF < C0ST3 ( IFACE , J, X) . EO. BEST1 ( IFACE) > THEN 
PRBCF(E) aPRBCP<K ) -1 
NOTES < J ) aNOTXS < J )- 1 
IF(PRBCP(X> . EQ. 0)THEN 

DO  3E0  I  lal , NOTAVL 
DO  310  .Mai, NOTH 
COST3  <11  ,  <M ,  X)  all  CM 
310  CONTINUE 

310  CONTINUE 

END  IF 

IF (NOTES < J) . EQ . 0 ) THEN 

DO  340  II-l.NOTAVL 
DO  330  KXal , NOPRB 
COSTS <  1 1 , J  ,  XX ) aB I GH 
3  30  CONTINUE 

340  CONTINUE 

END  IF 

VRITE(7,*)J. IFACE, X 
COTO  303 
ENDIF 

330  CONTINUE 
300  CONTINUE 
303  DO  300  J«l,NOTR 
DO  307  Kal , NOPRB 
COST3< IFACE. J,X)-SICN 
307  CONTINUE 
300  CONTINUE 
NN-NN-1 

I F ( NN . CT . 0 ) COTO  370 

REWIND  7 

HaO 

DO  410  lal, NOTRKS 
DO  400  Jal.NOTREQd) 

NaM«l 

DO  300  Xal.NOTAVL 
READ (7,*) ITS, IAR,IPRB 
IF (I AX. IQ. HITMEN 
C 

C  CONFUTE  TANKER  OFFLOAD  CAFABILITT 

C 

CALI  CIRCLE(CFLAT(H> ,CPL0M(H> , EARLAT(H) , EARLON(N) ,S,X) 
DISTNCal 

TIHSa (DISTNC/400 . ) *00 . 

FFa. 534«< (ARCF6V< ITB, IAR)«EARGVT< IAR, IFRI1 >/3 . )*07. 71 
FUELaTIHE*FF» . 001 

FUELOF (H ) aARCFCVC ITB, UR)-EAR6VT( IAR, IPRD-FUEL 
C 

C  CONFUTE  10HBBR  ONLOAD  CAFABILITT 

C 

CALL  CIRCLK BRLAT< I > , IILOMI I > , CFLAT(N) , CFLON(H) , I, T) 
ARDIST(H) al 
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IF ( J . EQ .  I) THEN 

TIHERai (  ARDIST(H) -103 . ) /TAS) *80 . 

NSBCaTXMSR/ 30. 

TLEFTaTIMER-NSEG*30  . 

SVa47«. 

DO  370  Lai, NSEC 

SVaCV- ( ( 10  . Z7+ . 7 1 8*CW) *  .  03 > 

370  CONTINUE 

BGVTC  PaCV- <(10 . 27+ . 7 1 S*GV) *TLEFT/ 1000. ) 

ELSE 

CALL  CIRCLE(£ARLAT<M-1) , EARLON(H-l) , CPLAT(H) ,CPLON(H> ,1,7) 
DISSTal 

TIHESaOISST/TAS'OO. 

NSSCaTXMES/  30. 

TLEFTaTIHES-NSEG*  30 . 

GVaOOO. 

DO  300  Lai, NSEC 

CVaCV- ( ( 10 . 27+ . 7 1 8*CV) *.03) 

310  CONTINUE 

1CWTCP-CW-<<10.27+.718»CW)«TIEFT/1000. ) 

END  IF 

rUELa(00.33+<BGVrCF+4tf. >/2. *.  7037) / 1000 . ‘TIME 
BONBON <N) -410 . -BGVTCP+FUEL 
REVIND  7 
GO  TO  400 
ENDir 

300  CONTINUE 

400  CONTINUE 
410  CONTINUE 

C 

C 

C  ARCPS  AND  EARS  ARE  NOV  ADJUSTED  A  CERTAIN  DISTANCE 
C  IN  AN  ATTEMPT  TO  EQUATE  BONIER  ONLOAD  CAPABILITT  WITS 
C  WITH  TANKER  OFrLOAO  CAPABILITT.  THE  DISTANCE  THE  ARCP 

C  IS  MOVED  DEPENDS  ON  HOV  MUCH  THESE  TVO  NUMBERS  DIFFER 
C  INITIALLY.  NEV  BOMBER  TIMES  TO  THE  ARCFS  ARE  ALSO 

C  CALCULATED. 

C 

C 

401  JaO 

DO  420  I  a 1 , NOTRKB 
DO  41S  Kal.NOTREQM) 

JaJ«l 

Ea FVE LOF < J )- BONBON < J ) 

IF(ASSd) . CT. . 4) THEN 
IFd.GT.  0)THEN 

IFIE.LT. . 5)THEN 
ARDISTt J)aARDIST< J)+l . 

ELSEIFd.LT.  1 .  )TNEN 
ARDISTC J)aARDIIT< J>+10 . 

ELSEIFd.LT.  2.  JTHEN 
ARDISTC J)aARDIST( J)> 13 . 
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ELSEI F  d . LT. 3 . )THEN 
ARDISTf J) .ARDISTt  J) +  30  . 

ELSEIFd  .  LT.  4.  )THEN 
ARDISTt J>. ARDISTt J) +40 . 

ELSEI F (X . LT. 7 . )THEN 
ARDISTt J) .ARDISTt J) +85 . 

ELSEIFd.  LT.  10  .  )  THEN 
ARDISTt  J ) -ARDISTt  J) +100 . 

ELSEIFd. LT.  13  .  >THEN 
ARDISTf J) .ARDISTf J)+Z00. 

ELSEIFd  .  LT.  30.  )THEN 
ARDIST< J) .ARDISTf J) +3 00 . 

ELSEIFf X . CT. 30 . )THEM 
ARDISTt J). ARDISTf J) +400 . 

END  IF 

ELSEI F(X . LT. 0 . )THEN 

IF(ABSd)  LT.  .  3 )  THEN 
ARDISTt J). ARDISTf J) -4. 

ELSEIF(ABSfX) .LT. 1. )THEN 
ARDISTf J>- ARDISTt J>-«. 

ELSEIFf ABSfX) .IT. 2.  ITHEM 
ARDISTf J). ARDISTf JJ-14 . 

ELSEIFf ABSfX) .LT. 3 . )THEN 
ARDISTf J)« ARDISTt J)-Z8  . 

ELSEIF  t ABSfX ) . LT. 4 . ITHEM 
ARDISTf J)«ARDlSTtJ>-38. 

ELSEIFf ABSfX) .LT. 7. )THEM 
ARDISTf J). ARDISTf J>-IZ. 

ELSEIFf ABSfX) . LT. 10 . >THEN 
ARDISTt J). ARDISTf J>-»3 . 

ELSEIFf ABSfX) . LT. 15 . )THEN 
ARDISTt J) .ARDISTf J) -115 . 

ELSEIFf ABSfX) .LT.30 . )THEN 
ARDISTf J) .ARDISTf J) -200. 

ELSEIFf ABSfX)  CT. 30. )THEN 
ARDISTf J). ARDISTf J)-SI0. 

ENDIF 

END  IF 
ENDIF 

ENDDI  StJ)>ARDISTt<))  +  130. 

IFfX.CT.DSNDDISfJ). ARDISTf  J)  +  I3. 

IFfK.CT. 2 >ENDDISt J) .ARDISTf J) +10. 

CALL  LATLONf BBLATfl ) , BBLONt I ) , ARDISTf J) , COURSE f I ) ,8 ,T) 
CFLATf J)-S 
CFLONf J).T 

CALL  LATLONf BBLATf I) .BBLONt I) , ENDDI S f J > , COURSE f I ) , S , T) 
EARLATf J ) .8 
EARLONf J).T 

FESTINt  J) .  t  f  ARDISTf  <3 )  -104  . )  *10  .  )  /TAB 
415  CONTINUE 
4X0  CONTINUE 


C  THE  REMAINDER  OF  THE  PROCRAM  RECOMPUTES  THE  TANKER 
C  OFFLOAD  CAPABILITY  AND  BOMBER  ONLOAD  CAPABILITY  AND 

C  REITERATES  THIS  PROCEDURE  UNTIL  THE  TWO  NUMBERS  ARE 

C  WITHIN  400  POUNDS  OF  EACH  OTHER.  BOMBER  ENTRY  POINT 

C  FUEL  IS  COMPUTED  FOR  EACH  BOMBER,  AS  WELL  AS  THE  TOTAL 

C  ENTRY  POINT  FUEL  FOR  THE  BOMBER  FLEET.  ALSO,  THE 

C  DIFFERENCES  BETWEEN  BOMBER  FUEL  REQUIRED  AND  ACTUAL 

C  BOMBER  FUEL  IS  CALCULATED  FOR  EACH  BOMBER. 

C 

C 

REWIND  7 
C 

Na  0 

DO  >10  I al , NOTRKS 
DO  9Z0  Jal ,NOTRIQ< I) 

H-M+l 

DO  130  Kal , NOTAVL 
READ<  7 , *  > ITS, 1 AR, I PRB 
IF ( IAR . EO . M)THEN 
C 

C  COMPUTE  TANKER  OFFLOAD  CAPABILITY 
C 

CALL  CIRCLE(TBLAT( ITB) ,TBLON< ITB) ,CPLAT(M) ,CPLON(H> ,X,Y> 
OISTNCal 

TIMER-< (DISTNC-184 .  >*80.  > /TAS 
NSECaTIMER/30. 

TLEFTaTIMER-NSEG* 30 . 

CW-TKCWT(ITB) 

DO  700  Lai, NSEC 

FFa.B30*GW+Z . S3 

IF (CW . LE . ISO . >FF-.813*CW+10 . 8 

FUEL-FF* . 03 

GWaGW-FUSl 

700  CONTINUE 

FF- . 83B*GW*Z . 83 
IFIGW.LE.IOO. )FF-.013*GV+10.S 
CWaCW- ( <  FF«TLEFT) * . 00 1 > 

C 

C  DETERMINE  IF  HOLDINC  IS  REQUIRED 
C 

HTIHEaFESTIN(H) -TIMER* . 001 
NSEG-KTIHE/30. 

TLEFTaHTIME-NSEC*  30 . 

DO  710  LLal.NSEC 

FP-.»033*CW-18.» 

CWaCW- FF*  .  81 
710  CONTINUE 

FF- . 80Z9*CW-Z0 . SO 
CWaCW-FF*TLEFT/1800 . 

ARCPCWC ITB,H> -CW 
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C  COMPUTE  NEV  EABGVT  FOR  ASSIGNED  PRB  AND  NEV  EAR 
C 

CALL  CIRCLE (EARL AT (M) , EARLON(H) , PRBLAT( IPRB) , 

(  PRB  LON  (  IFR8  )  »I .  Y) 

GV.FRBCVT 
IF (X. IE. 140 )THEN 
TLEFTaX/TAS*80 . 

GVaGV* . 157«TL£FT 

ELSE 

DISTXaX-73. 

TIME.DI STX/TAS*  6  0 . 

NSECaTIME/30 . 

TLEFTaTIME-NSEG* 30 . 

DO  720  KKal , NSEC 

CVa(GV+.3241>/.9756 
720  CONTINUE 

GV.GV*<<  .  61 3 * ( GV+ 1 . 3 ) ♦ 10 . 8  )  *TLEFT/  1000  .  >*2. 

END  IF 

EARGUT(M, IPRB) «CW 
C 
C 

CALL  CIRCL£( CPLAT(H) , CPLON(M) , EARLAT(N) , EARLON(N) ,X, Y> 
DISTNC-X 

TIME  a ( DI STNC/ 400 . >*60. 

FF-.524*<<ARCPGV<ITB, I AR > ♦EARGVTt IAR, IPRB) ill. >♦67.71 
FUEL«TINI*FF*. 001 

FUELOF(M) aARCPGVI ITB, !AR)-EARCVT( IAR, IPR3)-FUEL 
C 

C  COMPUTE  BOMBER  ONLOAD  CAPABILITY 
C 

CALL  CIRCLE <BBLAT< I) ,BBLON<I> ,CPLAT(N) , CPLON(M) , X, Y) 

ARDIST(M) -X 
IF(J.EQ.  DTHEN 

TIMERa ( ( ARDIST(N) -103 . ) IT AS) *  60 . 

NSKCaTIMER/30 . 

TLEFTaTIMER-NSEC* 30. 

GWa*76. 

DO  640  L-l ,NSEC 

GVaGV- ( <10.27+. 718*GV) * . 03  > 

640  CONTINUE 

SCVTCPaGV- < ( 10 . 27* . 7 1 8  *GV) *TLEFT/ 1000 . > 

ELSE 

CALL  CIRCLE ( EARLAT(H- 1 ) , EARLONCM- 1 > , CPLAT(N) ,CPLON(H) ,X,Y) 
DISBT.X 

TIMESaDIS8T/TAS*60 . 

NSEC. TIMES/30. 

TLEFTaTIMES-NSEG*  30 . 

GV.486 . 

DO  630  lal.NSSG 

GVaGV- ( < 10 . 27+ . 716*CV)* . 03) 

•30  CONTINUE 

BCVTCPaCV-<  < 10 . 27+ . 710*GV) *TLEFT/ 1000 . > 

END  IF 
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FUEL-(t0.33+( 8GVTCP+48S .  >/2. * . 7117 > / 1000 . *TINE 
BOMBON(N) -488 . -BCVTCP+FUEL 
REWIND  7 
CO  TO  820 
END  IF 

830  CONTINUE 
820  CONTINUE 
910  CONTINUE 
C 
C 
C 

DO  330  I-l.NOTAVL 

YY  <  I  > -ABS ( FUELOF ( I ) - BONBON ( I >  > 

330  CONTINUE 
YYY-YYU) 

DO  340  J.2.NOTAVL 
IF(YY(J)  CT. YYY) YYY-YY(J) 

340  CONTINUE 

IF (YYY.CT. . 4 ) GOTO  401 

FSVM-0 

J-0 

DO  380  I-1.NOTHK3 
J- J+NOTREQ( I ) 

IF (NOTREQC I ) . EO . 8 ) THEN 
EPFACT t I ) -UNBEF ( I > 

ELSE 

CALL  CIRCLE(EFLAT( 1 ) , EFLONf I ) , EARLATI J) , EARLONC  J) , E, Y) 
DISTT.X 

TIMER-DI 8TT/TAS*S0 . 

NSEC-TINER/ 30 . 

TL  EFT-TINER-NSEC  *  3 0 . 

GV-488 . 

DO  330  K.l.NSEG 

IFICV.CT. 340. >FF-10. 27+.718»CV 
IFICV.LE . 340 . ) FF -40 . 73+ . 829*CV 
CW-CW-FF* .03 
330  CONTINUE 

IFIGV.GT. 340 . )FF-10 . 27+ . 718*GV 
IF(CV.LE . 340 . >FF-4I . 73+ . 028»CV 
GW-CW-FF*TLEFT/ 1000 . 

EFFACTI 1 )«CV-21 0 . 3 
END  IF 

FSUH.FSUM+EFFACTd) 

EXDIF  < I ) .EPFACT < I ) -BFFREQ ( I > 

300  CONTINUE 
ITXN-ITSR+l 
!F< ITER. LE. 3) GOTO  ISO 
END 
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SUBROUTINE  Cl RCLE ( LAT, LONG , ELAT, ZLONG , 01 8T, COURSE) 

C  THIS  SUBROUTINE  COMPUTES  THE  GREAT  CIRCLE  COURSE  AND  DISTANCE 
C  BETWEEN  TWO  POINTS 

COMMON / FACTOR / P I ,  RAO 

REAL  LAT, LONG , 01  ST, COURSE , ELAT, ZLONG 

Pl-LAT/RAO 

PIMaLONC/ RAO 

P 2. EL AT /RAD 

PZM-E LONG /RAD 

00-1.570710327 

IF(PI.CT.QD)  Pl-QD 

IF<P2 .GT.QD)  PZaQO 

OaACOS(SIN( PI >  *SIN< PI ) *COS( PI ) *COS( PZ) *COS( P2H-P1M) > 

RHO-D*  3437. 748 77 

THETA-ACOS ( (SIN(PZ)-SIN(Pl) *COS( D) ) /SIN<D) / COS ( PI ) ) 

IF(SIN(PZM-P1M) .GEO)  THETA- 2*P I -THETA 

DIST-RHO 

COURSE -THETA* RAD 

RETURN 

END 

C 

C 

SUBROUTINE  LATLONt LAT, LONG , 01  ST, COURSE , ELAT, ZLONG ) 

C 

C  THIS  SUBROUTINE  COMPUTES  A  NEW  LATITUDE  AND  L0NC1TVDE  CIVEN 

C  AN  INITIAL  LATITUDE  AND  LONGITUDE  AND  A  DISTANCE  AND  COURSE. 

C 

COMMON/ FACTOR / P I , RAO 

REAL  LAT, LONG , DIST, COURSE , ELAT, ZLONG 

RKO-OIST 

Pl-LAT/RAO 

PIN-LONG /RAO 

THETA - COURSE /RAD 

RaRHO/3417. 74077 

P2-ASIN(SIN(P1) *COS( R ) +C0S (P1)*S1N(R) * COS (THETA) ) 

ZEE- ( COS (R ) -SIN(P1)*SIN(P2))/C0S(P1) /COS(P  2) 

D-ACOS(ZEE) 

IF(SIN(THETA)  GE.O.Q)  D--D 
P2H-P1H+D 

IF( P2H. CE . 0.0. AND . SIN< PZM) . LT. 0 . 0 >  P2H-P2H-2*PI 

IFIPZM.LT. 0 . 0 . AND.SIN(PZM) . GT.0.0)  PZM-PZM«Z*PI 

ELAT-PZ*RAD 

ELONCaPZM*RAD 

RETURN 

END 
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Introduction 


This  appendix  describes  some  of  the  characteris¬ 
tics  of  the  GNET  subroutine.  It  contains  excerpts  from 
References  2  and  7.  The  actual  subroutine  is  not  included 
in  this  report  because  it  is  proprietary  information.  The 
reader  is  directed  to  Reference  2  for  information  on  how 
to  obtain  the  actual  program  listing. 

The  Code  and  Its  Capabilities 

GNET  is  a  machine  independent  FORTRAN  program  for 
the  solution  of  the  capacitated  transshipment  problem. 

The  capacitated  transshipment  problem  is  the  most  general 
of  the  minimum  cost  flow  models  which  include  the  capaci¬ 
tated  and  uncapacitated  transportation  problems  and  the 
personnel  assignment  problems. 

The  capacitated  transshipment  model  and  its  special¬ 
izations  are  minimum  cost  network  flow  problems.  The  goal 
is  to  determine  how  (or  at  what  rate)  a  good  should  flow 
through  the  arcs  of  a  network  to  minimize  shipment  costs. 

The  transportation  and  assignment  models  are  simplifications 
of  this  transshipment  formulation. 

Unlike  the  usual  textbook  approach,  data  is  stored 
only  for  the  arcs  that  are  present  in  the  network.  This 
saves  storage  and  computations  since  for  most  practical 
problems  every  node  is  not  connected  to  every  other  node. 


It  is  also  permissible  to  use  multiple  arcs  to  model  piece- 
wise  linear  convex  shipping  costs. 

GNET  performs  a  primal  network  simplex  algorithm 
by  structural  manipulation  of  a  list  structure  representa¬ 
tion  of  the  network  minimization  problem  and  its  triangu¬ 
lated  bases.  Each  basis  of  the  network  problem  is  stored 
as  a  rooted  arborescence  and  the  pivotal  transformations 
are  performed  structurally  rather  than  by  numerical  matrix 
operations.  GNET  has  been  designed  with  selective  safe¬ 
guards  that  guarantee  successful  optimization  from  both 
programming  and  mathematical  viewpoints.  The  code  exhibits 
the  following  general  features: 

A.  Exact  Solutions.  The  solutions  produced  are 
absolutely  free  of  rounding  error.  (Data  and  all  calcula¬ 
tions  are  integer.) 

B.  Problem  Size.  The  code  has  been  calibrated, 

4 

tuned,  and  tested  on  problems  with  nominal  sizes  of  10 
nodes  and  105  arcs.  Modified  versions  have  successfully 
solved  problems  an  order  of  magnitude  larger. 

C.  Machine  Independence.  The  routine  is  written 
in  FORTRAN  V  and  has  been  tested  and  tuned  on  most  major 
computer  systems  including  AMDAHL,  Burroughs,  CDC,  Honey¬ 
well,  IBM,  UN I VAC,  and  TI/ASC.  The  user  must  specify 
the  largest  representable  absolute  integer  ("BIGI”)  and 
make  adjustments  only  for  any  nonstandard  FORTRAN 
restrictions. 


131 


D.  Storage  Economy.  The  algorithm  uses  much  less 
memory  for  execution  than  any  known  competing  method. 

Using  IBM  "words"  as  a  unit,  the  net  region  requirement 
for  (standard  version)  GNET  and  data  arrays  is  1,600  4-  9M 
4-  3N  where  M  is  the  number  of  nodes,  and  N  is  the  number 
of  arcs. 

E.  Non-Cycling  Algorithm.  The  algorithm  cannot 
cycle  in  the  presence  of  primal  degeneracy.  A  terminal 
solution  is  guaranteed. 

F.  Robustness.  GNET  has  been  tuned  on  hundreds 
of  diverse  problem  formulations.  The  candidate  queue 
pricing  mechanism  provides  an  extremely  powerful  problem- 
adaptive  control  for  pivot  trajectories  which  automatically 
exploits  special  problem  structure.  The  pricing  parameters 
for  the  candidate  queue  are  automatically  set  at  default 
values  for  excellent  general  performance.  However,  signifi¬ 
cant  further  improvements  in  efficiency  are  possible  by 
custom  adjustment  of  the  parameters  for  classes  of  problems 
exhibiting  very  unorthodox  structure. 

G.  Adaptability.  GNET  can  solve  network  problems 
of  extremely  large  size  without  further  modification. 
Moreover,  the  candidate  queue  mechanism  is  designed  to 
permit  modification  for  truly  huge  problems. 
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Subroutine  NTRD  Input  Example 

GNET  uses  an  input  subroutine  called  NTRD.  The 
subroutine  reads  in  the  input  data  from  a  file  using  the 
following  format  and  structure: 

1.  Header  card  (15)  number  of  nodes  in  the 
network,  M. 

2.  Arc  description  cards: 


column 

format 

item 

1-4 

A4 

"Name"  of  arc  (optional) 

5-6 

2X 

(Not  used) 

7-12 

16 

Source  of  node  of  arc 

13-18 

16 

Destination  node  of  arc 

19-20 

2X 

(Not  used) 

21-30 

no 

Cost  per  unit  flow 

31-40 

no 

Capacity  of  arc 

41-50 

no 

Lower  bound  of  arc 

End  of  file. 

Supplies  and  demands  can  be  input  by  using  "dummy" 
arcs  (not  explicitly  stored  after  input) .  NTRD  will  assume 
two  "artificial  nodes"  numbered  M+l  and  M+2.  Node  M+l 
is  a  "super  source”  to  the  problem,  and  node  M+2  is  a 
"super  sink."  Supplies  to  the  network  are  specified  as 
"dummy”  arc  capacities  from  node  M+l  to  the  source  nodes 
involved.  Demands  are  given  as  "dummy"  arc  capacities  from 
the  destination,  or  sink  nodes,  to  node  M+2. 
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All  costs  must  be  integers.  Use  o£  decimal  points 
on  input  is  not  allowed.  Also,  all  integers  must  be 
"right-oriented"  in  the  appropriate  fields.  Zero  costs 
are  admissible. 

All  arcs  must  have  capacities  that  are  nonnegative 
integers.  To  create  an  "uncapacitated"  arc,  give  the  arc 
a  positive  capacity  just  greater  than  the  total  of  all 
supplies  to  the  problem. 

Nodes  should  be  numbered  1,2,...,M  with  positive 
integers.  It  is  not  necessary  to  use  all  the  numbers 
between  1  and  M.  Numbering  is  arbitrary  with  respect  to 
arc  orientation.  Any  sequence  of  arc  description  records 
is  acceptable.  Multiple  arcs  between  any  two  nodes  are 
admissible. 

GNET  solves  the  all-equation  transshipment  problem, 
thus  total  supply  must  equal  total  demand.  Inequality 
problems  are  easily  transformed  to  this  form  by  the  addi¬ 
tion  of  an  extra  node  and  "slack"  arcs.  The  input  is 
described  here  for  the  usual  case  of  zero  lower  bounds, 
positive  capacities  and  no  arcs  with  fixed  flow. 
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Sample  Problem  Input 

COLUMNS :  1234 56 789^1234  56  78 9^1234 56  789^1234 56 789^1234 56 789§ 

header :  7 


arcs:  AB 

1 

3 

3 

100 

0 

PC 

1 

4 

6 

100 

0 

PC 

1 

2 

0 

50 

0 

BC 

3 

4 

1 

200 

0 

BD 

3 

7 

1 

300 

0 

CD 

4 

7 

2 

100 

0 

CF 

4 

5 

6 

150 

0 

I>- 

7 

9 

0 

120 

0 

FT) 

6 

7 

11 

200 

0 

E- 

6 

9 

0 

70 

0 

PC 

5 

4 

5 

150 

0 

FE 

5 

6 

8 

90 

0 

EE 

5 

6 

5 

50 

0 

GC 

2 

4 

4 

50 

0 

GF 

2 

5 

2 

400 

0 

-A 

8 

1 

0 

90 

0 

-G 

end:  *eor 

8 

2 

0 

100 

0 

Problem  Feasibility  and 
Scaling  of  Arc  Costs 

the 

For  some  minimum  cost  flow  models,  such  as  the 


classic  assignment  problem,  there  is  always  a  feasible 
solution  and  thus  am  optimal  solution.  In  general,  however 
even  if  total  supply  equals  total  demand,  there  may  not  be 
a  feasible  solution. 

GNET  introduces  artificial  arcs  into  the  initial 


basis  and  then  drives  their  flows  to  zero  for  a  feasible 
solution.  At  termination,  if  the  flow  on  all  artificial 
arcs  is  zero,  the  final  solution  is  feasible,  and  thus 


optimal.  If  any  artificial  arc  has  positive  flow,  the 
problem  has  no  feasible  solution. 

For  the  artificial  arcs,  GNET  calculates  a  large 
positive  cost,  BIGM,  which  guarantees  that  the  flows  on 
these  arcs  will  be  zero  for  a  feasible  solution.  However, 
the  size  of  the  problem,  the  scale  of  the  costs,  and  the 
value  of  the  largest  integer  representable  on  the  machine 
may  limit  BIGM  to  something  less  than  this  calculated 
value.  If  BIGM  is  limited  in  this  way  a  warning  is  printed 
and  a  cost  scaling  flag,  ISCALE,  is  set  equal  to  1.  In 
this  case,  GNET  will  still  terminate  normally:  if  the  final 
solution  is  feasible,  it  is  also  optimal;  if  the  final 
solution  is  not  feasible,  it  is  possible  (although  this  is 
highly  unlikely  because  the  criterion  is  so  conservative) 
that  the  problem  actually  does  have  an  undiscovered  feasible 
solution.  Rescaling  the  costs  to  a  smaller  range  will  avoid 
this  potential  problem. 
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The  survivability  of  the  strategic  bomber  force  during 
Emergency  War  Order  missions  is  of  primary  concern  to  the  Stra¬ 
tegic  Air  Command.  Since  the  Soviet  Union  and  its  Warsaw  Pact 
allies  possess  the  most  extensive  air  defense  networks  in  the 
world,  a  penetrating  bomber  force  must  fly  as  low  as  possible  for 
as  long  as  possible.  This  tactic  minimizes  probability  of 
detection  and  vulnerability  to  defensive  threats.  It  also  greatly 
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increases  the  fuel  required  to  complete  the  mission.  This  addi- 
tional  fuel  is  supplied  by  one  or  more  in-flight  refuelings. 

The  initial  objective  of  this  thesis  was  to  develop  a  method 
for  assigning  tankers  to  the  bomber  force  in  an  optimal  manner. 

As  the  study  progressed  however,  it  became  clear  that  obtaining 
a  truly  optimal  solution  using  mathematical  programming  techniques 
cannot  be  guaranteed  due  to  the  nature  and  complexity  of  the 
problem.  As  a  result  the  emphasis  of  the  study  was  shifted  to 
developing  an  improved  method  for  solving  the  problem. 

Two  heuristic  methods  were  investigated.  The  first  method 
used  network  theory  in  an  attempt  to  minimize  the  costs  of  assign¬ 
ing  tankers  to  the  bombers.  The  second  method  was  based  on  the 
so-called  "greedy"  method.  This  method  basically  made  the  assign¬ 
ments  in  the  order  of  decreasing  marginal  cost  improvements. 

These  two  methods  were  evaluated  against  each  other  and  the  cur¬ 
rent  method  by  means  of  several  example  problems.  Both  methods 
yielded  better  results  than  the  one  currently  in  use,  with  the 
network  method  appearing  to  be  the  best. 
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